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FOREWORD 


The  1995  version  of  the  aeroprediction  code  (AP95)  was  limited  to  zero  roll  orientation 
(^>  =  0  deg  for  fins  in  plus-fin  arrangement)  aerodynamic  calculations  of  weapons  with  axisynunetric 
bodies.  The  AP95  had  nonlinear  total  aerodynamic  loads  computed  to  angle-of-attack  (AO A)  of 
90  deg,  but  these  loads  were  not  distributed  over  the  wings  and  body  for  more  effective  use  by 
structural  analysts.  Moreover,  the  AP95  predictions  of  axial  force  at  high  AOA  also  needed 
improvement.  As  a  result  of  these  limitations,  new  technology  was  developed  over  the  past  3  years 
to  eliminate  these  deficiencies.  These  new  technologies  were  integrated  into  the  AP95  and  the  new 
version  of  the  code  will  be  transitioned  as  the  AP98. 
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1.0  INTRODUCTION 


1.1  REQUIREMENTS 

Estimating  missile  aerothermodynamics  over  the  flight  regime  where  missiles  fly  is  quite 
important  in  all  phases  of  design.  These  aerodynamics  are  used  by  the  flight  dynamicist  to  estimate 
range  performance  and  miss  distance;  the  heating  information  is  used  to  perform  heat  transfer 
analysis;  and  the  aerodynamic  and  thermodynamic  loads  are  used  by  the  structural  engineer  to 
estimate  structural  integrity  of  the  configuration.  Missiles  that  are  launched  from  a  vertical  launcher 
can  experience  Angles  of  Attack  (AOA)  approaching  90  deg  if  a  strong  crosswind  is  present. 
Missiles  that  are  launched  from  aircraft  undergoing  maneuvers  can  also  experience  AOA 
approaching  60  deg.  Finally,  terminally  guided  missiles  undergoing  maneuvers  in  the  endgame  can 
anticipate  AOA  as  high  as  40  deg.  Hence,  it  is  necessary  to  have  aerodynamics  estimates  that  cover 
a  broad  range  of  flight  conditions  in  terms  of  Mach  number  and  AOA.  In  addition  to  M„  and  AOA 
requirements,  aerodynamics  are  generally  needed  as  a  function  of  roll  orientation.  This  is 
particularly  true  for  a  six-degree-of-freedom  (6DOF)  flight  dynamics  model.  It  is  also  true  for  a 
three-degree-of-freedom  (3DOF)  trim  flight  dynamics  model.  A  trim  model  is  defined  as  that  AOA 
where  pitching  moment  is  zero.  Most  of  the  world’s  missiles  fly  in  a  roll-stable  position  of 
^>  =  0  deg  (fins  in  “+”  fin  arrangement)  or  ^>  =  45  deg  (fins  in  “x”  fin  arrangement).  As  AOA 
increases  above  about  10  to  20  deg,  trim  aerodynamics  of  the  missile  configurations  in  these  two 
roll-stabilized  positions  can  be  quite  different. 

In  addition  to  the  freestream  requirements  of  AOA,  roll,  and  Mach  number,  several  geometry 
requirements  are  also  necessary.  The  first  of  these  is  associated  with  the  requirements  for 
axisymmetric  bodies.  Here,  the  nose  can  be  sharp,  blunt,  or  truncated  with  various  shaped  ogives 
present.  The  afterbody  is  typically  cylindrical  of  various  lengths  and  may  or  may  not  have  a  boattail 
or  flare  present.  Figure  1  is  a  sketch  of  a  body-alone  axisymmetric  case.  One  or  two  sets  of  lifting 
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surfaces  can  also  be  present  on  the  configuration  as  shown  in  Figure  1 .  The  cross-section  of  these 
lifting  surfaces  is  typically  symmetric  with  a  sharp  or  blunt  biconvex,  double-wedge  or  modified- 
double-wedge  airfoil  section. 


The  desire  to  increase  weapon  range  and  maneuverability,  to  design  weapons  that  are  more 
optimum  from  an  aircraft  total  drag  and  radar  signature  standpoint,  or  to  provide  optimum  loadout 
of  multiple  missiles  in  a  ship’s  vertical  launcher  has  driven  weapons  designers  to  consider 
nonaxisymmetric  body  shapes.  Some  typical  shapes  are  shown  in  Figure  2.  While  most  missiles  in 
the  United  States  and  foreign  countries  in  existence  today  have  axisymmetric  body  configurations, 
these  conceptual  design  tradeoffs  of  various  configurations  other  than  axisymmetric  require 


FIGURE  2.  NONCIRCULAR  CROSS-SECTION,  WING-BODY  CONFIGURATIONS  OF 
INTEREST  FOR  THE  AEROPREDICTION  CODE 
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engineering  estimates  of  aerodynamics.  Current  state-of-the-art  methods  for  predicting 
aerodynamics  of  nonaxisymmetric  body  shapes  with  engineering  accuracy  are  much  more  limited 
than  for  axisymmetric  bodies.  This  is  primarily  driven  by  the  fact  that  to  get  reasonable  accuracy 
of  the  aerodynamics  requires  an  accurate  description  of  the  body  geometry.  To  describe  the 
geometry  of  a  complex  body  shape  accurately  can  take  days  or  weeks  depending  on  the  requirements 
of  the  aerodynamics  code  being  used. 


1 .2  LIMITATIONS  OF  1995  VERSION  OF  THE  AEROPREDICTION  CODE  (AP95) 

The  latest  version  of  the  aeroprediction  code,  AP95‘,  has  several  limitations  when  viewed 
against  the  overall  requirements.  First  of  all,  the  code  is  limited  to  the  roll  position  of  <^  =  0  deg 
aerodynamics.  Since  die  code  is  a  semiempirical  model  used  for  preliminary  design,  it  is  important 
that  aerodynamics  be  available  in  the  roll  position  of  O  =  45  deg  as  well.  This  new  technology  was 
therefore  developed^  and  will  be  a  part  of  the  AP98.  Secondly,  while  the  AP95  defined  nonlinear 
aerodynamic  loads  to  AOA  90  deg,  it  did  not  distribute  all  these  loads  over  the  body  and  lifting 
surfaces  to  make  the  code  more  useful  to  the  structural  engineer  who  might  be  using  the  code  to 
perform  a  beam  analysis  of  a  missile  structure.  This  problem  was  also  addressed  and  documented 
in  Reference  3. 

A  third  problem  that  existed  with  the  AP95  was  the  accuracy  of  the  axial  force  coefficient 
at  high  AOA.  This  problem  became  more  prominent  at  subsonic  Mach  numbers,  where  the  axial 
force  could  go  negative,  and  at  all  Mach  numbers  when  a  control  deflection  was  of  opposite  sign  to 
the  AOA.  New  technology  was  developed'^  to  correct  this  problem  and  this  new  technology  will  also 
be  a  part  of  the  AP98.  The  final  weakness  of  the  AP95  compared  to  the  aeroprediction  requirements 
is  the  limitation  of  axisymmetric  bodies.  This  was  an  extremely  difficult  problem  to  resolve. 
However,  after  a  couple  of  years  where  several  difficult  technology  shortfalls  were  overcome,  this 
limitation  has  also  been  removed  and  documentation  on  the  new  technology  has  just  been 
completed.^ 

Figure  3  is  a  summary  that  shows  the  natural  evolution  of  the  aeroprediction  code  from  the 
first  version  in  1972  to  the  seventh  and  most  recent  version,  the  AP98.  This  figure  shows  how  major 
new  technology  was  developed  to  meet  the  tactical  weapon’s  requirements.  The  remainder  of  this 
report  will  summarize  the  theoretical  methods  contained  in  the  aeroprediction  code;  compare  the 
aerodynamic  predictions  of  the  latest  version  to  those  of  the  AP95  and  to  experimental  data;  and, 
finally,  summarize  and  briefly  state  what  new  technology  is  still  needed  in  the  realm  of  approximate 
aerodynamic  codes  for  future  development. 


2.0  LOW  ANGLE  OF  ATTACK  (AOA)  AERODYNAMIC  METHODS 


Most  of  the  low  AOA  aerodynamic  methods  are  state-of-the-art  and  have  been  used  for  many 
years.  They  are  primarily  theoretical  in  nature  and  based  on  small  perturbation  theory  or  local  slope 
or  hypersonic  flow  methods.  They  have  been  well  documented  in  past  references.  As  a  result,  each 
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topic  will  be  briefly  summarized  with  appropriate  references  given  for  those  who  want  to  pursue  the 
methods  more  thoroughly. 


FLIGHT  CONDITIONS 


VERSION 

WEAPONS 

AERODYNAMIC.S 

MACH 

REAL  GAS 

AOA  RANGE 

ROLL 

NONLINEAR  COMPUTERS 
DISTRIBUTED 

LOADS 
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AVAILABLE 

1972 

AXISYMMETRIC 

UNGUIDED 

PROJECTILES 

STATIC  ONLY 

0-3 

NO 

o 

O 

o 

O 

ii 

NO 

CDC 

1974 

AXISYMMETRIC 

MISSILES 

PROJECTILES 

ROCKETS 

SAME 

0-3 

NO 

SAME 

SAME 

NO 

CDC 

1977 

SAME 

STATIC  AND 
DYNAMIC 

0-3 

NO 

SAME 

SAME 

NO 

CDC,  IBM 

1981 

SAME 

SAME 

0-8 

NO 

0-15“ 
(LIMITED 
CONF.  AT 
HIGHER  a) 

SAME 

NO 

CDC,  IBM,  VAX 

1993 

SAME 

SAME 

0-20 
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0-30“ 

SAME 

NO 

CDC,  IBM,  VAX, 

SILICON 

GRAPHICS 

1995 

SAME 

SAME 

SAME 

YES 

o 

1 

VO 

o 

SAME 

NO 

INTERACTIVE 

PC 

1998 

AXISYMMETRIC 

AND 

ASYMMETRIC 

SAME 

SAME 

YES 

SAME 

^  =  0“,  45“ 

YES 

IMPROVED 

INTERACnVE 

PC 

MISSILES, 
PROJECTILES, 
AND  ROCKETS 


FIGURE  3.  EVOLUTION  OF  AEROPREDICTTON  CODE  IN  TERMS  OF  MAJOR  NEW  ADDED  CAPABIUTY 


2. 1  BODY  ALONE  METHODOLOGY 


21.1  Hybrid  Theory  of  Van  Dyke  Combined  with  Modified  Newtonian  Theory 

This  method  was  first  introduced  in  Reference  6  and  was  the  key  to  being  able  to  obtain 
accurate  pressure  distributions  and  axial  and  normal  force  coefficients  and  center  of  pressure  on 
blunt,  truncated  or  sharp  nosed  bodies  of  revolution.  Reference  6  was  thus  the  foundation  of  the  first 
version  of  the  aeroprediction  code  (AP72)  and  to  the  author’s  knowledge,  remains  to  this  day  the 
only  accurate  way  (short  of  a  numerical  method  to  solve  the  full  equations  of  motion)  to  obtain  static 
aerodynamics  in  the  Mach  number  range  of  1.2  to  1.8.  Reference  6  combined  the  Hybrid  Theory  of 
Van  Dyke’’^  (HTVD)  with  the  MNT.®  The  key  to  the  accurate  computation  of  pressure  coefficients 
all  along  the  surface  of  a  blunt  nosed  body  was  the  starting  solution  and  match  point  between  the 
MNT  and  HTVD. 
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The  HTVD^  combines  a  second-order  solution  to  the  potential  equation  with  a  first-order 
crossflow  solution  first  espoused  by  Tsien.*  The  advantage  of  this  method  is  that  it  gives  second- 
order  accuracy  in  the  axial  direction  where  first-order  accuracy  is  generally  unacceptable  for  drag 
computations.  On  the  other  hand,  first-order  accuracy  in  the  crossflow  plane  is  typically  acceptable 
for  normal  force  and  center  of  pressure  computations.  The  fundamental  reason  for  this  is  that 
perturbations  in  the  flow,  due  to  the  presence  of  a  body,  have  more  impact  in  the  axial,  as  opposed 
to  the  normal  force  direction.  Hence,  to  get  axial  force  accuracy  compatible  with  a  goal  of 
±10  percent  requires  second-order  methods;  whereas  ±10  percent  accuracy  on  can  be  obtained 
with  first-order  methods,  in  many  cases. 

The  Hybrid  theory  comes  from  the  potential  equation  of  fluid  mechanics.  It  is  limited  to 
supersonic  flow  (we  have  used  this  method  down  to  M«  =  1.2)  where  the  assumption  of  isentropic 
flow  (shock  waves  are  weak)  can  be  made.  This  t5q)ically  limits  the  upper  Mach  number  range  to 
about  M„  =  2.0  to  3.0,  depending  on  the  body  shape.  Also,  the  slope  of  the  body  surface  must  be  less 
than  the  Mach  angle.  The  Tsien  solution,  or  crossflow  part  of  the  solution,  comes  from  the 
linearized  perturbation  equation.  On  the  other  hand,  the  second-order  solution  to  the  axial  flow  is 
found  by  obtaining  a  particular  solution  to  a  reduced  version  of  the  full  potential  equation.  This  is 
the  key  to  the  accuracy  improvement  afforded  by  Van  Dyke’s  solution  in  that  some  of  the 
nonlinearity  inherent  in  the  axial  flow  problem  is  brought  into  the  solution  by  this  process.  The 
beauty  of  the  Van  Dyke  method  is  that  this  particular  second-order  solution  is  given  entirely  in  terms 
of  the  first-order  solution.  That  is,  one  simply  solves  the  first-order  perturbation  solution  for  the 
axial  flow  and  then  solves  an  algebraic  equation  for  the  second-order  solution  where  the  boundary 
condition  at  the  body  is  satisfied. 

Newtonian  Impact  Theory  assumes  that,  in  the  limit  of  high  Mach  number,  the  shock  lies  on 
the  body.  This  means  that  the  disturbed  flow  field  lies  in  an  infinitely-thin  layer  between  the  shock 
and  body.  Appl5dng  the  laws  of  conservation  of  mass  and  momentum  across  the  shock  5delds  the 
result  that  density  behind  the  shock  approaches  infinite  values  and  the  ratio  of  specific  heats 
approaches  unity.  The  pressure  coefficient  on  the  surface  becomes® 

Cp  =  2sin26^  (1) 


where  6^^  is  the  angle  between  the  velocity  vector  and  a  tangent  to  the  body  at  the  point  in  question 
(see  Figures  2  and  3).  is  defined  by: 

sin  (6gq)  =  sin0  cosa  -  sina  cost})  cos0  (2) 

As  noted  in  the  earlier  discussion,  the  Hybrid  Theory  is  limited  to  conditions  where  the  body 
slope  is  less  than  the  local  Mach  angle.  This  means  it  is  not  applicable  in  the  nose  region  of  a  blunt 
missile.  On  the  other  hand,  MNT  gives  very  acceptable  estimates  of  pressure  coefficients  in  the  nose 
region,  even  for  low  supersonic  Mach  numbers  where  the  assumptions,  inherent  in  the  Newtonian 
Impact  Theory,  are  violated.  As  already  mentioned,  the  key  to  the  successful  combination  was  in 
the  starting  solution.  At  low  supersonic  Mach  numbers,  the  pressure  overexpands  on  a  blunt  nose 
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tip  as  it  proceeds  around  the  blunt  portion  from  the  stagnation  point  to  a  given  portion  of  the  nose. 
In  order  to  capture  this  overexpansion.  Reference  6  found  that  it  was  necessary  to  start  the  HTVD 
near  its  maximum  acceptable  slope  and  allow  the  pressure  to  expand  around  the  surface. 
Simultaneously,  the  MNT  was  started  at  the  stagnation  point  and  allowed  to  expand  until  the 
pressure  coefficients  of  the  MNT  and  the  HTVD  were  equal.  This  was  defined  as  the  match  point. 
Upstream  of  the  match  point,  MNT  was  used  in  the  force  and  moment  calculations,  whereas 
downstream,  HTVD  was  used. 


2.1.2  Second-Order-Shock-Expansion  Theory  (SOSETl  Combined  with  MNT 

First-order  Expansion  Theory  was  first  proposed  by  Eggers  et  al.  for  bodies  of  revolution 
flying  at  high  supersonic  speeds.^®  Basically,  the  Shock-expansion  Theory  computes  the  flow 
parameters  at  the  leading  edge  of  a  two-dimensional  (2-D)  surface  with  the  oblique  shock  wave 
relations  and  with  the  solution  for  a  cone  at  the  tip  of  a  three-dimensional  (3-D)  body.  Standard 
Prandtl-Meyer  Expansion  (PME)  is  then  applied  along  the  surface  behind  the  leading  edge  or  tip 
solution  to  get  the  complete  pressure  distribution  over  the  body  surface.  This  theory  inherently 
assumes  that  the  expansion  waves  created  by  the  change  in  curvature  around  the  body  are  entirely 
absorbed  by  the  shock  and  do  not  reflect  back  to  the  body  surface.  Since  the  theory  assumes  constant 
pressure  along  one  of  the  conical  tangent  elements  of  the  surface,  fairly  slender  surfaces  must  be 
assumed  or  many  points  along  the  surface  assumed  to  obtain  a  fairly  accurate  pressure  distribution. 
Another  way  of  stating  this  is,  to  minimize  the  strength  of  the  disturbance  created  by  Mach  waves 
emanating  from  the  expansion  comer  and  intersecting  the  shock,  the  degree  of  turn  should  be  small. 

Syvertson  et  al.  extended  the  generalized  Shock-expansion  Theory  on  pointed  bodies  and 
sharp  airfoils  to  what  he  called  a  second-order  theory.”  He  defined  the  pressure  along  a  conical 
frustum  by 


P  =  Pc  -  (Pc  -  P2)  e’’*' 


(3) 


instead  of  a  constant  on  each  segment  as  was  the  case  in  the  generalized  theory.  Here,  p^  is  the 
pressure  on  a  cone  with  the  given  cone  half  angle  equal  to  the  slope  of  the  conical  segment  with 
respect  to  the  axis  of  sjmrimetry.  The  term  p2  is  the  pressure  just  aft  of  a  conical  segment  which  is 
calculated  from  a  PME  of  the  flow  around  a  comer. 

Also,  Til  is  a  decay  constant  that  must  be  positive.  If  it  is  negative,  it  is  set  to  zero  and  the 
second-order  theory  reverts  back  to  the  first-order  theory. 

Jackson  et  al.”  combined  SOSET  with  MNT  to  treat  blunt-nosed  configurations  with  or 
without  flares.  Jackson  et  al.”,  like  Syvertson  and  Dennis,”  assumed  that  the  lifting  properties  could 
be  predicted  by  assuming  that  the  original  body  is  made  up  of  several  equivalent  bodies  of  revolution 
represented  by  the  various  meridians.  They  assumed  the  match  point  between  the  MNT  and  second- 
order  shock  pressure  prediction  to  be  the  angle  that  corresponds  to  shock  detachment  on  a  wedge 
with  the  given  freestream  Mach  number. 
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DeJamette  et  made  significant  improvements  to  the  work  of  Jackson  et  al.^^  and 
Syvertson."  These  new  improvements  included  the  following: 

1.  An  exact  (as  opposed  to  an  approximate)  expression  for  the  pressure  gradient 
downstream  of  a  comer. 

2.  A  new  expression  for  pointed-cone  pressures  at  angle  of  attack  that  improves  the  initial 
pressure  prediction  over  that  of  the  tangent  cone  theory. 

3.  A  new  technique  for  calculating  pressures  of  bodies  at  incidence. 

The  pressure  computations  at  AO  A,  showed  improvement  over  the  method  of  Jackson.*^ 
DeJamette,  et  al.*^  derived  a  new  expression  for  pointed-cone  pressure  at  a  >  0  by  combining  Slender 
Body  Theory  (SBT),  Newtonian  Theory,  and  an  approximate  expression  for  Cp  to  give: 

a=0 

Cp  (a,  0,  <J),  M)  =  Cp^^^  +  ACp  (4A) 


where 


ACp  =  -sin2a  sin20  coscj) 


sin^a  cos^0 


2-ij  (l-tan20)-^2+|-j  sin^t}) 


(4B) 


Cp  =  sin^0. 


,  (Y  ^  1)  .  2  p  J  Y  ^  1  .  J_' 

(Y  -  1)  +  2  [  2 


(4C) 


and 


=  (M„^  -  1)  sin20^ 


Note  also,  that  while  Equation  (4)  was  strictly  defined  for  pointed  cone  pressures  at  angle  of  attack, 
it  could  also  be  used  in  a  Tangent  Cone  sense  to  obtain  pressures  at  any  point  on  a  body  surface. 
DeJamette  actually  used  loading  functions  to  obtain  body  alone  lift  properties,  however.*^ 

Equations  (1)  through  (4)  are  thus  the  methods  used  to  calculate  pressures  between  the.  Mach 
numbers  of  about  1 .8  and  6.0.  Below  M.  of  1 .8,  HTVD  combined  with  MNT  is  used.  Above  M. 
of  6.0,  a  slight  improvement  of  Equation  (4B)  was  fovmd  to  be 
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ACp  =  -  (2a)  sin  (20)  cos  (4))  +  (Fcos^0)  +  (4/3  sin  (20)  cos  ((|))  )  a^  (5A) 


where 


F  =  (2  -  -i)  (1  -  tan^  0^)  -  (2  +  sin2(J) 


or 


,  _  (2a)  sin  (20)  cos  (6) 

,p  =  -  - - ^ -  (53) 


Equation  (5  A)  is  used  for  pointed  body  configurations,  as  well  as  for  blunt  body  configurations  in 
the  windward  plane  area  60  deg  <  (j)  £  180  deg.  For  the  leeward  plane  area  on  blunt  bodies, 
Equation  (5A)  is  replaced  by  Equation  (5B). 


2.1.3  Skin  Friction  Drag  Methodology 

Four  options  are  allowed  in  the  AP98  for  skin-friction  drag  computations.  These  are:  (1)  all 
laminar  flow,  (2)  wind  tunnel  model  with  no  boundary  layer  trip,  (3)  wind  tunnel  model  with 
boundary-layer  trip  (all  turbulent  flow),  and  (4)  typical  flight  conditions.  The  code  assumes 
transition  from  laminar  to  turbulent  flow  of  1  x  10®  and  0.5  x  10®  for  the  body  and  fins  respectively 
for  a  typical  flight  vehicle  that  has  about  normal  roughness.  For  a  wind  tunnel  model  with  no 
boundary  layer  trip,  these  values  of  transition  Reynolds  numbers  have  been  increased  to  4  x  10®  and 
2  X  10®  respectively.  Also,  the  values  are  changed  (decreased)  as  angle-of-attack  increases  to  allow 
turbulent  flow  to  occur  on  a  larger  portion  of  the  vehicle. 

For  the  portion  of  the  vehicle  where  turbulent  flow  occurs,  the  Van  Driest  method  is  used 
for  skin-friction  drag  computations.  Reference  1  gave  a  method  for  computing  the  laminar  portion 
of  the  skin-friction  drag.  The  reader  is  referred  to  References  1,15,  and  16  for  more  details  of  the 
skin-friction  drag  methodology  used  in  the  AP98. 


2.1.4  Base  Drag 

The  AP72®  developed  a  simplified  method  to  compute  base  drag  on  bodies  of  revolution 
using  an  average  of  several  data  bases  of  base  pressure  coefficient  on  cylindrical  afterbody 
configurations  as  a  function  of  freestream  Mach  number.  Turbulent  flow  ahead  of  the  base  was 
assumed.  Boattail  effects  were  accounted  for  as  well  as  a  first  approximation  to  AOA  effects  in 
Reference  6.  The  AP74^^  added  fins  to  the  body  so  as  to  be  able  to  compute  aerodynamics  on 
missiles  and  rockets  as  well  as  unguided  projectiles.  As  a  result,  a  rough  approximation  to  the  body 
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base  drag  as  a  result  of  tail  fins  was  made.  Also,  the  fins  themselves,  if  they  had  blunt  trailing  edges, 
had  a  base  drag  of  their  own.  This  was  based  on  a  two-dimensional  base  pressure  coefficient  data 
base. 


Later  on,  a  request  was  made  to  the  National  Aeronautics  and  Space  Administration/Langley 
Research  Center  (NASA/LRC)  to  provide  additional  base  pressure  data  to  allow  a  more  thorough 
estimate  of  the  effects  of  AO  A,  fin  thickness  and  location  effects,  as  well  as  Mach  number.  These 
tests  were  conducted  by  NASA/LRC;  an  improved  empirical  method  for  predicting  base  drag  was 
developed  and  documentation  of  the  tests  and  methodology  was  given  in  Reference  17. 

Figure  4  summarizes  the  empirical  base  drag  prediction  method  of  Reference  17.  The 
parameters  Fj,  F2,  and  F3  of  Figure  4  are  each  defined  by  a  set  of  curves  in  Reference  17.  Once  A 
or  B  of  Figure  4  has  been  obtained,  then  power-on  or  boattail  (or  both)  effects  are  computed  in  an 
additive  manner.  While  the  model  of  Reference  17  has  been  developed  for  power-off  and  power-on 
effects'*  computed  independently,  the  methodology  of  combining  the  two  together  in  an  additive 
manner  is  an  assumption  and  has  never  been  adequately  validated.  It  is  suspected  that  the  Figure  4 
methodology  could  be  improved  upon  with  the  combined  effects  of  power-on  in  conjunction  with 
fins  and  AO  A  effects,  if  this  type  of  data  were  available. 


2.1.5  Real  Gas  and  Aeroheating  Effects 

The  main  reason  the  fourth  version*^  of  the  aeroprediction  code  was  limited  to  Mach  number 
8  was  that,  above  M„  =  6,  real  gas  effects  start  becoming  important  but  can  still  be  neglected  at 
M„  =  8.  However,  as  Mach  number  increases  substantially  above  M^  =  6,  the  need  to  include  real 
gas  effects  into  the  aeroprediction  code  increases  if  one  is  interested  in  inviscid  surface  temperatures. 
If  one  is  only  interested  in  forces  and  moments,  real  gas  effects  have  a  slight  effect  on  the  pitching 
moment,  but  only  second-order  effects  on  axial  and  normal  force.^  However,  one  of  the  key  issues 
in  high-speed  vehicles  is  aerodynamic  heating,  material  selection,  and  insulation.  Any  excess  weight 
can  have  a  strong  adverse  impact  on  vehicle  performance.  Thus,  a  simple  yet  accurate  method  of 
estimating  vehicle  surface  temperature  (inviscid)  for  use  in  heat  transfer  analysis  is  needed. 

To  extend  the  Aeroprediction  Code  to  Real  Gases  required  the  SOSET  to  be  extended  to 
include  real  gas  solutions.  This  required  a  real  gas  solution  for  a  cone;  a  Prandtl  Meyer  expansion 
for  real  gases;  a  derivation  of  a  pressure  derivative  (which  was  part  of  the  “n,  in  Equation  (3);  and  a 
way  to  compute  inviscid  temperature,  knowing  pressure.  This  was  accomplished  and  documented 
in  Reference  14. 

The  inviscid  surface  temperature  was  not  the  main  interest  area  in  extending  the 
aeroprediction  code  to  real  gases.  The  main  interest  was  in  providing  better  estimates  of  heat  transfer 
rate;  heat  transfer  coefficients;  and  recovery  or  adiabatic  wall  temperature  at  each  computational 
point  along  the  vehicle  surface.  This  information  could  then  be  used  in  a  heat  transfer  code  as  the 
convective  part  of  the  heat  transfer.  This  is  typically  the  most  difficult  part  of  the  heat  transfer  to 
calculate.  Standard  SOTA  methods^'  were  used  in  obtaining  the  convective  portion  of  heat  transfer. 
These  methods  are  Fay-RiddelP^  on  the  blunt  nose,  Beckwith  and  Gallagher^*  on  the  blunt  leading 
edge  of  a  wing  and  the  Eckert  Reference  Enthalpy  method^''  on  the  remaining  body  and  wing  points. 
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Body  Alone 


F,  =  Body  Alone  a  Effects 


Body  with  Tail  Fins 

1 .  Deflection  and  Thickness  Effects 


F2=  |a  +  6|  Effects 

F3  =  Additional  Effects  due  to  Thickness 
2.  Fin  Location  Effects  when  x/c  *  0 

(Cp )  =  (Cp  \  +  0.01  (ACp  \ 

'  ®/a,  6,  t/c,  x/c  \  a  ' 


6,  t/c,  x/c 


Power-On 

ACp  due  to  Power-On  by  Modified  BrazzeP^  Method  Added  to  Above  Values  of 


Boattail 


C.,  -  -  Cp. 

Flare 

Ca,  =  -Cp/dp/d,,/ 


FIGURE  4.  EMPIRICAL  BASE  DRAG  PREDICTION  MODEL  OF  AP98 
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Also,  a  mass  balance  technique  was  used  to  account  for  the  fact  that  at  high  Mach  number,  the 
boundary  layer  swallows  the  entropy  layer.  As  a  result,  the  entropy  at  the  outer  edge  of  the  boundary 
layer  is  no  longer  the  value  at  the  surface  of  the  body,  but  an  adjusted  value. 


2.2  WING  ALONE  METHODOLOGY 


2.2.1  Lifting  Surface  Theorv^^ 

Lifting  Surface  Theory  refers  to  the  solution  of  the  flow  over  a  3-D  wing  where  the 
distribution  of  pressure  is  allowed  to  vary  in  both  the  spanwise  and  chordwise  direction.  The 
fundamental  equation  is  the  3-D  perturbation  equation,  along  with  boundary  conditions  that  require 
the  flow  to  be  tangent  to  the  wing  surface  and  the  velocity  on  the  upper  and  lower  surfaces  to  be 
equal  at  the  trailing  edge. 

The  assumptions  involved  in  the  Lifting  Surface  Theory,  as  applied  to  most  missile 
configurations,  are  therefore  small  perturbations  in  the  flow  due  to  the  presence  of  the  wing  and  the 
thickness  and  camber  effects  are  zero  or  small  compared  to  AOA  effects.  The  current  version  of  the 
aeroprediction  code  neglects  thickness  effects  in  calculating  the  lift  on  the  wings.  Also,  the  Prandtl- 
Glauert  rule^®  is  used  to  relate  the  compressible  subsonic  normal  force  and  pitching  moment  of  the 
wing  to  that  calculated  at  Mach  number  zero.  Using  the  Prandtl-Glauert  rule  allows  one  to  solve  La 
Place’s  equation  at  Mach  number  zero  and  then  compute  compressible  subsonic  wing-alone 
aerodynamics  at  any  Mach  number  less  than  about  0.8  to  0.9,  where  transonic  effects  become 
important. 

To  compute  roll  and  pitch-damping  derivatives,  the  flow  tangency  boundary  condition  is 
changed  to  relate  the  local  AOA  to  roll  or  pitch  rate.  A  more  detailed  discussion  of  the  wing-alone 
linear  aerod)mamics  is  given  in  References  27  and  28. 


2.2.2  Three-Dimensional  Thin  Wing  Theory  TTPYWri^^ 

TDTWT  is  quite  similar  to  lifting  surface  theory  (LST)  in  the  sense  the  same  perturbation 
equation  is  used.  The  only  difference  is  that  TDTWT  is  normally  used  to  represent  the  supersonic 
flow  solutions  of  the  perturbation  equation,  versus  LST,  representing  the  subsonic  solutions.  Since, 
for  supersonic  flow,  solutions  to  the  perturbation  equation  are  hyperbolic  versus  elliptic  for  the 
subsonic  case,  they  generally  are  easier  to  obtain.  This  is  because  no  upstream  influence  is  felt  by 
a  disturbance  at  a  given  point  on  the  wing  surface.  In  contrast,  the  subsonic  solutions  required  a 
matrix  inversion  at  each  wing  element  to  determine  the  unknown  coefficients  used  to  determine  the 
pressure  differential  from  lower  to  upper  surfaces.  On  the  other  hand,  the  assumptions  of  TDTWT 
are  the  same  as  for  LST.  They  both  assume  small  perturbations  in  an  isentropic  flow.  The  isentropic 
flow  assumption  means  no  shock  waves  are  allowed. 
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In  contrast  to  the  body  solutions  generated  by  Van  Dyke,  adequate  wing  solutions  can  be 
obtained  at  higher  Mach  numbers.  This  is  because  of  the  low  slopes  present  on  most  wing  planforms 
(thickness  is  generally  very  small),  the  wing  frontal  area  is  generally  less  than  10  percent  of  the  body 
frontal  area,  and  in  the  region  of  leading  edge  bluntness,  where  perturbation  theory  is  invalid. 
Modified  Newtonian  Theory  is  used  for  wave  drag  calculation.  While  thickness  is  neglected  for 
normal  force  and  pitching  moment  calculations,  it  is  not  neglected  for  axial  force  computations.  A 
similar  process  to  subsonic  LST  is  used  to  compute  the  dynamic  derivatives  of  roll  and  pitch 
damping  where  the  local  AO  A  is  assumed  to  be  a  function  of  either  roll  or  pitching  motion. 

For  transonic  flow,  which  here  is  considered  to  be  Mach  numbers  in  the  range  of  0.8  to  1.2, 
interpolation  using  the  TDTWT  and  LST  methods  is  done  using  a  process  similar  to  that  in  the 
aircraft  DATCOM.2® 


2.3  INTERFERENCE  AERODYNAMICS 

There  are  two  main  types  of  interference  dealt  with  in  the  low  AOA  methods.  The  first  has 
to  do  with  the  interference  effect  the  wing  has  on  the  body  or  the  body  has  on  the  wing.  The  second 
has  to  do  with  the  downwash  effect  of  a  tail  surface  in  the  wake  of  a  forward  lifting  surface.  Each 
of  these  interference  effects  will  be  dealt  with  separately.  Wing  to  wing  or  shock  wave  interference 
will  not  be  discussed  at  present,  but  will  be  alluded  to  in  the  discussion  on  nonlinear  aerodynamic 
methods. 


2.3.1  Wing-Body  and  Body-Wing  Interference  Effects 

To  better  understand  the  interference  lift  components,  it  is  instructive  to  examine  the  total 
normal  force  of  a  configuration  as  defined  by  Pitts  et  al.^°  This  is  given  by 

“  ^Ng  [(^W(B)  ^B(W))  ^  ^W(B)  ^(W))  ^w] 


The  first  term  in  Equation  (6)  is  the  normal  force  of  the  body  alone  including  the  linear  and 
nonlinear  components;  the  second  term  is  the  contribution  of  the  wing  (or  canard),  including 
interference  effects  and  control  deflection;  the  third  term  is  the  contribution  of  the  tail,  including 
interference  effects  and  control  deflection;  and  the  last  terms  are  the  negative  downwash  effect  on 
the  tail  or  body  due  to  wing-shed  or  body-shed  vortices.  The  K’s  represent  the  interference  of  the 
configuration  with  respect  to  AOA,  and  the  k’s  represent  the  interference  with  respect  to  control 
deflection.  Each  of  these  interference  factors  is  estimated^^  by  slender  body  or  linear  theory.^®  As 
such,  they  are  independent  of  AOA.  Figure  5  plots  the  SBT  values  of  Kb(w),  kw(B)  and  kB(W)- 


12 


NSWCDD/TR-98/1 


0.0  0.2  0.4  0.6  0.8  1-0 

r/s 

HGURE  5.  SLENDER  BODY  THEORY  INTERFERENCE  LEFT  FACTORS 


As  the  Mach  number  increases  supersonically,  SET  gives  values  of  Kb(w)  which  are  too  high 
if  the  wing  is  near  the  missile  rear.  This  is  because  much  of  the  carryover  lift  onto  the  body  is 
actually  lost  to  the  wake  of  the  vehicle.  Linear  theory  formulations  are  available  for  the  infinite  and 
no-afterbody  cases  to  replace  Figure  5  if  the  parameter 

PAR  (1  +  X)  [If  (mp)  +  1]  >  4  (7) 

Moore^’  then  linearly  interpolated  between  the  infinite  and  no-afterbody  cases  as  a  function  of  the 
area  covered  by  the  Mach  lines  to  obtain  Kb(w)  for  the  short  afterbody  case. 

Strictly  speaking,  the  methodology  discussed  here  is  limited  to  slender  bodies  with  triangular 
planforms  of  low  aspect  ratio.  Experience  has  shown  that  if  the  correct  value  of  wing-alone  lift  is 
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computed,  the  interference  factors  can  give  very  reasonable  results  for  wings  which  do  not  have 
triangular  planforms  or  even  have  low  aspect  ratio.  Reference  27  showed  how  an  engineering 
estimate  to  interference  lift  could  be  obtained,  even  for  planforms  such  as  that  shown  in  Figure  6A. 
The  actual  SET  configuration  is  that  shown  in  Figure  6B.  Since  most  of  the  interference  lift  occurs 
near  the  wing-body  juncture,  Reference  27  used  approximations  given  by  Equation  (8), 

[^w(B)ln  “  1  ([Kw(b)^i 


^B{W)ln  ([^W(B)^I  f^W(B)3l)^  (8) 


to  estimate  the  interference  factors  of  the  wing  in  Figure  6A.  G  in  Equation  (8)  is  the  ratio  of  the 
root  chord  of  the  wing  for  which  the  interference  factor  is  desired  to  that  of  the  wing  that  slender 
body  theory  assumes.  That  is 


'  ("r)l 


2.3.2  Wing-Tail  hiterference  Effects 

The  last  two  terms  of  Equation  (6)  are  also  interference  terms.  Cj^  is  the  lift  on  the  tail 
caused  by  the  vortices  shed  by  the  wing  or  canard  upstream.  is  tBe  negative  lift  on  the 

afterbody  due  to  wing-shed  vortices.  These  terms  are  also  calculated*mal5rtically  and  are  given  by: 

2lt  {ARV  (4  -  O  A„, 


-4r 

AwV. 


(10) 
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HGURE  6A.  WING  FOR  WHICH  INTERFERENCE  LEFT  IS  DESIRED 


KNOWN: 

Kw(B) 

Kb(W) 

I  > 

-ll 

HGURE  6B.  ASSUMED  SLENDER  BODY  REPRESENTATION 
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Here  i  is  the  tail  interference  factor  given  by  Pitts  et  al.,^°  and  T  is  the  strength  of  the  wing-shed 
vortex. 


The  APS  1  and  prior  versions  of  the  aeroprediction  code  (APC)®’^’  compute  based  on 

Equation  (9),  but  neglect  of  Equation  (10).  As  a  result,  some  overprediction  of  nonnal  force 
occurs,  and  generally  speaking!  the  vehicle  is  too  stable  as  well.  Later  versions  of  the  APC,  which 
include  nonlinear  effects,  inherently  account  for  this  term  indirectly  through  the  direct  use  of  wind 
tunnel  data  bases  to  estimate  wing-tail  interference  nonlinearities.  These  nonlinear  methods  will  be 
discussed  later  in  this  report. 


2.4  EMPIRICAL  METHODS^’’^**’ 

It  is  fair  to  wonder  why  approximate  linearized  aeroprediction  codes  are  defined  as 
semiempirical  with  all  the  theoretical  methods  discussed  so  far.  The  truth  is  that  while  these 
methods  allow  the  individual  component  forces  and  moments  to  be  calculated  fairly  rigorously  at 
a  given  Mach  number  or  AOA,  there  are  still  many  conditions  where  the  analytical  methods 
presented  previously  are  either  not  applicable  or  the  difficulty  in  applying  them  is  not  worth  the 
effort.  In  those  cases,  empirical  methods  are  generally  used.  The  combination  of  theoretical  and 
empirical  techniques  in  a  code  is  thus  why  they  are  called  semiempirical  codes.  A  few  examples 
where  empirical  methods  are  used  are  transonic  aerodynamics,  body-alone  subsonic  aerodynamics, 
rotating  band  or  protuberance  aerodynamics,  and  base  drag  of  the  body  and  lifting  surfaces  (which 
was  discussed  already).  There  are  actually  analytical  methods  available  for  transonic  aerodynamic 
computations.  However,  most  of  the  methods  are  inconsistent  from  a  computational  standpoint  with 
the  approximate  codes.  What  is  done,  in  many  cases,  is  to  use  the  sophisticated  analytical  tools^’-^  *® 
to  estimate  the  transonic  aerodynamics,  as  a  function  of  key  geometric  parameters,  then  to  include 
these  into  an  engineering  code  in  a  table  lookup  fashion.  Obviously,  for  a  vehicle  that  spends  a  large 
portion  of  its  time  in  the  transonic  flow  region,  0.8  <  M„  <  1 .2,  it  would  be  justifiable  to  use  a  more 
sophisticated  estimation  process. 


3.0  fflGH  AOA  AERODYNAMIC  METHODS 


To  calculate  aerodynamics  at  high  AOA  requires  either  a  sophisticated  computational  fluid 
dynamics  (CFD)  code^*’^^  or  a  semiempirical  APC  which  bases  the  nonlinear  aerodynamics  on  wind 
tunnel  data  bases.  Both  approaches  (CFD  and  semiempirical)  are  needed  in  the  design  process.  The 
semiempirical  methods  tend  to  be  less  accurate  than  the  CFD  methods,  but  are  much  less  costly  and 
much  faster  to  get  answers.  As  a  result,  they  tend  to  be  used  in  the  conceptual  and  preliminary 
design  stages  of  development.  On  the  other  hand,  the  CFD  approaches,  while  being  more  costly  and 
time  consuming,  are  also  more  accurate,  particularly  when  a  CFD  code  has  been  benchmarked 
against  a  set  of  existing  experimental  data.  As  a  result,  CFD  codes  tend  to  be  used  more  in  the  later 
phases  of  design,  to  supplement  wind  tunnel  data,  or  to  answer  a  technical  question  as  to  what  is 
occurring  in  the  overall  flow  field. 
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While  generic  wind  tunnel  data  bases  are  not  necessary  for  CFD  codes,  they  are  essential  for 
semiempirical  codes.  The  general  approach  of  a  semiempirical  code  is  to  try  to  back  out  the 
nonlinear  aerodynamic  terms  based  on  wind  tunnel  data  and  an  estimate  of  the  linear  aerodynamics. 
Secondly,  an  attempt  is  made  to  mathematically  model  the  nonlinearity  based  on  approximations  to 
the  experimental  nonlinear  aerodynamic  terms.  Thirdly,  this  model  is  compared  to  other  data  sets 
and  the  mathematical  model  refined. 

Much  of  the  nonlinear  aerodynamic  methodology  for  the  AP98  has  been  developed  recently; 
hence,  a  more  thorough  summary  of  this  methodology  will  be  given  in  comparison  to  the  low  AOA 
aerodynamic  methods  discussed  in  Section  2.0  of  this  report.  Section  3.0  will  deal  with  the  new 
methods  for  axisymmetric  bodies,  and  Section  4.0  will  cover  the  noncircular  body  cross-section 
methodology. 


3. 1  AXIAL  FORCE  COEFHCIENT 

The  new  methodology  for  computing  the  nonlinear  term  of  axial  force  is  documented  in 
Reference  33.  Basically,  Reference  33  assumed 


C 


A 


(11) 


where  represents  the  zero  AOA  axial  force  coefficient,  and  ,  the  change  in  axial  force  as 
AOA  increases  or  decreases.  is  estimated  quite  well  using  the  low  AOA  methods  discussed  in 
Section  2.0.  However,  the  methodology  for  treating  in  the  AP95*  could  be  refined  considerably, 

particularly  at  subsonic  Mach  numbers  where  can  acUially  go  negative.  The  other  area  where 
of  the  AP95  required  improvement  was  when  a  and  8  were  of  opposite  sign. 

The  new  methodology  of  Reference  33  treated  C  *  as  a  fourth-order  equation  in  AOA.  That 

IS 


C.  =  f(M,a)  =  Aa  +  Ba^  +  Ca^  + 


(12) 


The  constants  A,  B,  C,  and  D  were  evaluated  using  several  wind  tunnel  data  bases  for  body-alone 
and  body-tail  configurations.  For  wing-body-tail  cases,  the  body-tail  parameters  were  adjusted 
somewhat  based  on  wind  tunnel  data  on  wing-body-tail  cases  in  conjunction  with  the  AP95.  The 
constants  A,  B,  C,  and  D  are  defined  by: 
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A=  f(M,0) 

B=  -3.509  f'(M,0)  +11.005  f(M,30)  -2.757  f(M,60)  +0.41  f(M,90) 

C=  3.675  f(M,0)  -17.591  f(M,30)  +7.041  f(M,60)  -1.179  f(M,90) 

D=  -1.181  f(M,0)  +6.771  f(M,30)  -3.381  f(M, 60)  +0.752  f(M,90) 

where  f'(M,0),  f(M,30),  f(M,60)  and  f(M,90)  are  the  values  of  the  function  f(M,a)  or  its  derivative 
at  AOA  0, 30, 60  and  90  deg,  respectively.  These  values  for  body-alone,  body-tail,  and  wing-body- 
tail  cases  are  shown  in  Figure  7. 

To  model  the  change  in  axial  force  with  control  deflection  and  AOA,  the  following  equations 
are  utilized: 


(14A) 


sinS^  f(M,aT) 


(14B) 


where 


«w 

=  a 

+  6^ 

CC.J. 

=  a 

(15) 


It  was  found  that  when  a  and  6  were  of  the  same  sign,  f(M,a)  was  close  to  1 .0.  However,  when  a 
and  6  were  of  opposite  signs,  the  axial  force  contributions  due  to  control  deflection  were  too  high. 
As  a  result,  empirical  expressions  were  derived  for  f(M,a)  at  the  two  roll  positions  of  ^>  =  0  and 
$  =  45  deg.  These  results  are  given  in  Figure  8.  Figure  8  was  derived  by  use  of  the  AP98  in 
conjunction  with  wind  tunnel  data  where  the  controls  were  deflected.  Note  that  the  only  roll 
dependence  of  comes  from  control  deflection.  IF  there  are  no  control  surfaces  deflected,  then 
is  approximately  the  same  within  the  present  prediction  accuracy  for  $  =  0  deg  and  O  =  45  deg.  ° 


3.2  NORMAL  FORCE  COEFHCIENT 

3.2.1  Configuration  Aerodynamics 

There  are  basically  two  approaches,  to  the  author’s  knowledge,  for  including  configuration 
nonlinear  aerodynamics  into  the  normal  force  Equation  (6),  repeated  here  for  convenience. 
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A.  BODY  ALONE 


B.  BODY  TAIL 


C.  WING-BODY-TAIL 

HGURE?.  AXIAL  FORCE  AOA  VARIATION  PARAMETERS 
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[(^W(B)  *  ^B(W))  “  f^W(B)  ^B(W))  ^w] 

+  [(Kj(B)  ^B(T))  ^  (^•(B)  ^B(T))  '*'  ^Nb(v)  (6) 

The  first  method  was  the  Equivalent  AO  A  (EAOA)  method.  In  basic  form,  the  EAOA  method 
defines  an  AOA  where  the  wing-alone  normal  force  is  twice  the  single  fin  value.  That  is,  following 
Reference  34,  the  normal  force  on  a  single  fin  at  =  0  deg  roll  and  low  AOA  is 

^Np,B,  =  t^W(B)«  +  K(B)^  +  (^“)w]  (17) 


Reference  34  then  defined  an  EAOA  as 

“eq  ^W(B)“  ^W(B)^  ■*"  (18) 

The  term  (Aa)^  of  Equations  (17)  and  (18)  is  the  average  AOA  induced  on  the  fins  by  vortices 
(either  body-  or  wing-generated).  The  EAOA  was  later  extended^^’^*  to  include  local  Mach  number 
and  other  configuration  effects.  The  modeP  itself  treats  Kw(b)  as  a  nonlinear  term,  but  assumes  Kb(w) 
is  based  on  SET.  The  major  advantage  of  this  approach  is  that  it  allows  computation  of  individual 
fin  aerodynamics  at  any  roll  orientation  and  thus  can  be  used  to  predict  6DOF  aerodynamics.  That 
is,  not  only  C^,  C^,  and  X^p,  but  Cy,  C„,  and  Q  as  well.  It  has  been  shown  to  be  effective  in 
predicting  missile  aerodynamics  with  reasonable  accuracy  for  many  cases.  The  major  weak  points 
are  that  the  method  is  still  an  empirical  approximation  and  the  many  analytical  approximations  and 
assumptions  apparently  degrade  the  accuracy  in  predicting  aerodynamics  over  the  direct  use  of  data. 
Secondly,  the  concept  is  confusing  to  many  people.  That  is,  by  lumping  terms  together  as  in 
Equation  (18),  it  is  hard  to  follow  the  physics  of  each  aerodynamic  phenomena. 

As  a  result  of  these  two  shortcomings,  the  present  author  attempted  a  different  approach  to 
incorporating  configuration  nonlinear  aerodynamics  into  Equation  (6),  beginning  in  1992,”  and 
further  defined  in  References  1  and  2.  The  approach  of  References  37,  1,  and  2  was  to  take  each 
term  of  Equation  (6)  and  break  it  down  into  a  linear  component,  which  is  estimated  by  the  theoretical 
low  AOA  methods  of  Section  1,  and  a  nonlinear  term,  which  is  estimated  directly  from  several  large 
missile  component  data  bases.^*"^  Since  the  approach  at  NSWCDD  for  estimating  6DOF 
aerodynamics  was  either  numerical  codes^*’”  or  wind  tunnel  data,  it  was  decided  to  focus  on  trim 
aerodynamic  models  for  the  APC.  These  models  are  used  to  predict  planar  aerodynamics  of  weapons 
designed  to  fly  in  a  roll-stabilized  mode  of  either  0  deg  or  45  deg.  This  makes  the  job  of  getting 
accurate  aerodynamics  much  easier  than  predicting  6DOF  aerodynamics.  It  also  ties  the  physics  of 
the  flow  much  closer  to  the  mathematics  since  each  term  of  Equation  (6)  must  have  a  linear  and 
nonlinear  estimate.  Of  course,  the  disadvantages  of  the  direct  approach  compared  to  the  EAOA 
approach  are:  a)  one  must  derive  nonlinear  expressions  for  all  aerodynamic  terms  of  Equation  (6) 
versus  lumping  several  together  as  done  by  Equation  (18);  and  b)  it  has  not  been  derived  for  6DOF 
aerodynamics.  Using  the  approach  of  References  1  and  2,  we  have  been  able  to  obtain  average 
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accuracy  levels  of  ±10  percent  on  and  and  ±4  percent  of  body  length  on  Xcp.  By  “average” 
is  meant  enough  AOAs  or  Mach  number  aeropredictions  so  as  to  get  a  good  statistical  sample. 
Usually,  the  first  few  degrees  AOA  are  excluded  from  the  sample  due  to  measurement  errors. 

The  way  Equation  (6)  is  implemented  in  the  AP93,  AP95,  and  AP98  codes  is  by  computing 
the  nonlineaiities  of  the  body-alone  term  based  on  AOA  only,  but  all  the  other  nonlinear  corrections 
to  the  interference  factors  and  wing  or  tail  aerodynamics  are  based  on  the  total  AOA,  or  a  +  6.  This 
was  done  for  ease  of  implementation  into  the  operational  AP81  code.  Hence,  even  though  the 
component  nonlinear  terms  were  derived  from  wind  tunnel  data  bases  from  AOA  data  only,  they 
were  implemented  in  a  total  local  AOA  sense.  The  empirical  nonlinear  corrections  were  then  fine- 
tuned  based  on  comparisons  to  other  configuration  aerodynamics  outside  the  data  bases.  The  center 
of  pressure  of  individual  missile  component  aerodynamics  were  treated  similarly  to  the  AP81  except 
shifts  in  center  of  pressure  based  on  data  were  derived.  These  shifts  were  implemented  in  tabular 
form  as  a  function  primarily  of  AOA  and  M„. 

Also,  the  Cj^  term  of  Equation  (6),  which  is  the  downwash  normal  force  on  the  body  due 
to  the  wing-shed  vortices,  is  neglected.  This  is  because  it  is  inherently  included  in  the  wind  tunnel 
data  bases,  and  it  is  believed  the  errors  in  trying  to  analytically  estimate  the  term,  subtract  it  out  on 
one  configuration,  and  then  add  it  back  in  later  on  a  different  configuration,  are  as  large  or  larger  than 
the  errors  from  incorporating  it  into  the  term. 

Equation  (6)  can  also  be  rewritten  as 


C  =  C 


+  c 


N, 


'W(B) 


+  c 


N, 


'B(W) 


'N. 


T(B) 


+  c 


-Bcn 


N. 


T(V) 


(19) 


where  it  is  understood  that  Cj^  encompasses  the  term.  For  ease  of  implementation  into  an 
existing  code  designed  prim^S’y  for  linear  aerodynamics,  most  of  the  terms  in  Equation  (19)  are 
separated  into  a  linear  and  nonlinear  contribution  due  to  a  or  6.  For  example,  the  wing-body  term 
is  computed  in  the  AP98  code  as  follows: 


^W(B) 


(20) 


The  linear  or  small  AOA  terms  of  Equation  (20)  are  estimated  by  linear  theory  (LT)  or  SBT  (see 
Section  2).  This  gives  the  Aeroprediction  Code  a  good  fundamental  basis  for  its  aerodynamic 
estimates.  The  nonlinear  corrections  due  to  higher  AOA  or  control  deflection  are  each  estimated 
directly  from  component  wind  tunnel  data  bases.^*"*®  Each  of  the  other  terms  in  Equation  (19)  is 
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treated  in  a  similar  fashion  to  Equation  (20)  in  the  actual  implementation  into  the  Aeroprediction 
Code. 


In  the  context  of  Equation  (19),  we  therefore  seek  the  nonlinear  definition  of  each  of  the 
terms  in  Equation  (19).  It  is  expected  that  the  body-alone  term  [first  term  of  Equation  (19)]  will  be 
independent  of  (j).  In  reality,  this  is  not  necessarily  the  case  for  M  <  2  and  high  AOA  because  of  the 
asymmetric  shedding  of  vortices.  The  mechanism  of  this  shedding  is  not  clear,  but  it  is  suspected 
that  slight  imperfections  in  the  flow  or  body  shape,  from  uniform  or  axisymmetric,  respectively, 
could  contribute  to  this  phenomenon.  At  present,  the  Aeroprediction  Code  does  not  account  for  out 
of  plane  aerodynamics,  and  therefore  the  side  force  created  by  the  asymmetric  shedding  of  body 
vortices  is  not  predicted.  Also,  in  the  Reference  38  data,  normal  force  varied  by  about  10  percent 
as  a  function  of  roll  in  the  region  of  asymmetric  vortex  shedding.  Instead  of  including  this  variation, 
it  was  averaged  out. 

Since  the  AP98  includes  aerodynamics  at  both  the  =  0  and  ^  =  45  deg  roll  orientations, 
roll  dependence  on  wing-body  aerodynamics  were  investigated  based  on  slender  body  and  linear 
theory  A  detailed  summary  of  these  findings  is  presented  in  Reference  2.  A  summary  of  the  key 
findings  in  Reference  2  on  roll-dependent  aerodynamics  is  repeated  here  for  convenience. 

a)  For  cruciform  wings  alone  or  a  wing-body  combination,  the  total  normal  force  is 
independent  of  roll. 

b)  For  a  planar  wing-body  combination  at  roll,  the  loading  on  the  windward  plane  panel  is 
greater  by  an  equal  amount  to  that  on  the  leeward  plane  panel.  This  means  that  if  one 
were  trying  to  design  a  code  for  lateral  aerodynamics,  roll  dependence  of  each  fin 
planform  must  be  considered.  On  the  other  hand,  if  longitudinal  aerodynamics  are  of 
primary  interest,  the  total  normal  force  on  the  entire  wing  planform  can  be  considered. 

c)  For  a  craciform  wing-body-tail  configuration  at  roll,  eight  vortices  are  shed  in  the  wing- 
body  region,  which  adversely  affects  the  tail  lift.  This  is  as  opposed  to  four  vortices  at 

=  0  deg. 

d)  The  planar  theory  developed  for  wing-tail  interference  can  be  used  to  approximate  the 
loss  of  lift  on  the  tails  at  =  45  deg. 

e)  The  aerodynamics  of  a  cruciform  wing-body-tail  combination  with  zero  control 
deflections  are  independent  of  roll  position. 

These  findings  for  roll  dependence  from  linearized  or  slender  body  theory  are  quite  useful 
in  helping  plan  how  to  develop  a  nonlinear  APC  for  both  ^>  =  0  and  $  =  45  deg.  While  the 
conclusions  of  linear  theory  roll  dependence  may  not  translate  to  the  nonlinear  case,  we  will  still  use 
the  findings  to  help  guide  the  nonlinear  code  development.  In  particular,  the  item  (a)  conclusion 
implies  use  of  the  ^>  =  0  deg,  wing-alone  data  for  <I>  =  45  deg.  This  is  quite  important  because  the 
available  wing-alone  data  bases  are  all  at  3>  =  0  deg.  This  means  that  any  nonlinear  wing-alone  roll 
dependence  will  be  included  in  the  interference  factors  rather  than  the  wing-alone  solution,  which 
is  independent  of  4>. 
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The  second  major  result  of  the  key  slender  body/linear  theory  roll  dependence  findings  is  that 
for  cruciform  missiles,  we  can  use  the  same  interference  approaches  as  in  the  AP95,  except  the 
constants  need  to  be  changed  because  of  a  different  roll  angle.  The  combination  of  these  two 
conclusions  are  quite  important  because  they  basically  allow  the  direct  usage  of  the  AP95  code  with 
different  constants  for  the  nonlinear  interference  terms  at  <E>  =  45  deg  versus  $  =  0  deg. 

The  third  significant  conclusion  is  that  for  small  AOA,  wing-body-tail  aerodynamics  are 
independent  of  roll  position.  This  allows  the  usage  of  wing-tail  interference  methodology  designed 
for  planar  computations  for  different  roll  orientations,  so  long  as  the  proper  number  of  vortices  are 
considered.  Again,  different  nonlinear  corrections  are  expected  for  the  <E>  =  45  deg  versus  the 
0  =  0  deg  roll  position. 

With  the  approach  of  how  the  APC  incorporates  nonlinear  configuration  aerodynamics 
having  been  discussed,  the  next  few  sections  will  take  each  of  the  terms  in  Equations  (6)  and  (19) 
and  discuss  how  the  nonlinear  term  is  defined.  Section  2  of  this  report  has  already  discussed  how 
the  low  AOA  aerodynanndcs  were  calculated. 


3.2.2  Body  Alone 

The  nonlinear  component  of  the  body-alone  normal  force  is  predicted  by  the  Allen-Perkins^^ 
viscous  crossflow  theory.  In  developing  this  theory,  Allen  reasoned  that  the  total  force  on  an 
inclined  body  of  revolution  is  equal  to  the  potential  term  discussed  in  Section  2,  plus  a  crossflow 
term.  This  term  is  based  on  the  drag  force  experienced  by  an  element  of  a  circular  cylinder  of  the 
same  diameter  in  a  stream  moving  at  the  cross  component  of  the  stream  velocity,  V„  sin  a.  This 
crossflow  term  is  primarily  created  by  the  viscous  effects  of  the  fluid  as  it  flows  around  the  body, 
often  separating  and  creating  a  nonlinear  normal  force  coefficient.  In  equation  form,  the  so-called 
viscous  crossflow  theory  is: 


i 


\ 


sin^a 


(21) 


Here,  ti  is  the  drag  proportionality  factor  or  crossflow  drag  of  a  cylinder  of  finite  length  to  one  of 
infinite  length.  is  the  crossflow  drag  coefficient.  has  two  possible  values  in  the  crossflow 
Mach  number  range  from  about  0. 1  to  about  0.6  depending  on  whether  the  flow  is  supercritical  or 
subcritical.  The  criteria  for  determining  if  the  flow  is  supercritical  is  that 


(22) 


where 
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fcosacota  ^  ^  \  J 

1  +[1  +(‘/2Cota)2]’'^' 

\  *  ^  olii  lA  /  1 

2  1  2  11 

[1  + 

Mv,  =  M„  sin  a 

M  00 

The  term  Rj^  in  Equation  (22)  is  Reynolds  number  based  on  body  diameter.  The  standard  values 
of  Rj^  and  l5lc  set  in  the  AP98  are  180,000  and  0.1  respectively.  However,  the  user  is  allowed  to 
change  these  values  if  a  known  value  of  Rj^^  and  Me  are  available  from  a  given  wind  tunnel  test. 

The  AP98  begins  transition  from  subcritical  to  supercritical  flow  when  the  Equation  (22)  is 
satisfied.  However,  the  minimum  value  of  C^j  in  Equation  (16)  is  not  obtained  until 

Rf^  =  Rn  +  25,000 

This  way,  the  flow  transitions  from  subcritical  to  supercritical  flow  in  a  continuous,  versus  abrupt, 
way. 


The  final  change  to  the  AP95  methodology  is  to  take  into  account  the  fact  that  bodies  with 
fins  present  have  different  physical  flow  characteristics  than  bodies  alone.  Most  of  the  wind  tunnel 
data  upon  which  the  characteristics  of  Rj^  and  M^^  are  determined  is  based  on  body-alone  data. 
Since  in  the  vicinity  where  wings  or  tails  are  present,  the  mechanism  of  boundary  layer  reattachment 
in  the  leeward  plane  will  be  harder  to  establish,  and  as  a  result,  the  length  of  the  body  where  the 
lifting  surfaces  are  present  will  be  taken  out  of  the  area  for  the  minimum  value  of  .  This  length 
is  assumed  to  be  two  root  chord  lengths.  Thus,  if  (C^  )j  and  (C^  ^  are  sub-  and  supercritical  values 
of  Cj  respectively,  then  a  modified  nonlinear  normm  force  of  the  body  alone  for  supercritical  flow 
(where  wings  are  present)  is: 


2c j  ^ 

The  total  normal  force  of  the  body  alone  is  then  the  sum  of  the  linear  term  from  Section  2  and  the 
nonlinear  term  from  Equation  (23). 


2c. 


(23) 


3.2.3  Wing  Alone 

One  of  the  primary  reasons  for  analyzing  the  slender  body  and  linear  theory  implications  on 
component  aerodynamics  in  Reference  2  was  to  show  that  the  wing-alone  methodology  developed 
for  the  #  =  0  deg  plane  could  also  be  applied  in  the  O  =  45  deg  plane.  This  was  at  least  trae  in  the 
linear  sense.  It  will  also  be  assumed  to  be  true  in  the  nonlinear  sense.  Any  nonlinearities  not 
accounted  for  in  using  the  =  0  deg  methods  for  =  45  deg  will  therefore  be  included  in  the 
interference  factors. 
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As  a  result  of  this  approach,  the  methodology  of  Reference  16  can  be  used  directly.  This 
methodology  uses  a  fourth-order  equation  in  AOA  to  estimate  wing-alone  lift,  as  opposed  to  the 
second-order  approach  of  Reference  1 1.  The  specific  equations  are  defined  by 

+  a2“w  +  *  aX  (24) 

a,  =  34.044  -  4.824  *  0.426  (c4_^.  -  6.412a,  <35) 

a,  =  -88.240  .  23.032  - 

>4  =  53-219  W..,,.  -  '■'■595  * 


^•522  (C,)„^.  4  U.464a,  ^26) 


2.661  -  5.971a,  (27) 


The  term  aj  of  Equations  (24)  through  (27)  is  the  value  of  wing-alone  lift  curve  slope  at  cc  =  0  given 
by  linear  theory.  The  terms  of  the  wing-alone  normal 

force  coefficients  at  oc  =  15, 35,  and  60  deg,  respectively,  defined  by  the  data  bases  of  References  39, 
40,  and  43.  Above  of  60  deg,  extrapolation  of  the  aerodynamics  at  a„„  of  60  deg  is  used.  For 
more  details  of  the  method,  the  reader  is  referred  to  Reference  44.  As  shown  in  Reference  44,  the 
fourth-order  method  of  Equations  (24)  through  (27)  improves  the  wing-alone  prediction  accuracy 
below  a,,  =  30  deg  over  the  second-order  method  used  in  the  AP93,  while  allowing  wing-alone 
aerodynamic  estimation  to  180  deg. 

The  value  of  Cj^  computed  by  Equation  (24)  includes  the  linear  and  nonlinear  term.  To 
include  it  into  a  term  like  liquation  (20),  requires  this  term  to  be  separated  into  a  linear  and  nonlinear 
component.  This  is  easily  done,  as  the  linear  term  is  known;  so  the  nonlinear  term  is  simply  the 
difference  between  the  total  wing-alone  value  and  its  linear  counterpart.  The  secant  slope  is  then 
formed  for  each  of  these  terms  by  dividing  the  linear  and  nonlinear,  parts  by  the  local  AOA,  a  +  6. 
This  approach  was  taken  for  ease  of  incorporation  into  an  existing,  operational,  linearized  code. 


3.2.4  hiterference  Aerodynamics 

Interference  aerodynamics  treated  in  a  nonlinear  sense  includes  that  associated  with  K^^), 

Equation  (6).  The  linear  form  of  these  interference  effects  are  all 
estimated  based  on  slender  ^dy  or  linear  theory  and  were  discussed  in  Section  2  earlier.  Basically, 
the  first  four  interference  factors  are  all  estimated  using  the  general  form 

K  =  Ksbt  +  AK(M^,  a,  AR,  A) 
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The  nonlinear  term  AK  of  Equation  (28),  is  defined  in  terms  of  five  tables  that  are  functions  of  Mach 
number,  AOA,  aspect  ratio  and  taper  ratio  for  both  Kw(B)  and  Kb(w).  The  general  nonlinear  trend  of 
those  two  interference  terms  is  shown  in  Figure  9.  This  general  trend  is  basically  the  same  for  both 
the  ^  =  0  and  =  45  deg  roll  orientations.  However,  the  five  tables  of  data  are  different  for  $  =  0 
and  $  =  45  deg  for  both  Kw(B)  and  Kb(W).  As  already  discussed,  for  cruciform  missiles,  SET  gives 
roll  independence  for  the  low  AOA  vaJues  of  Kw(b)  and  Kb(W).  Values  of  the  tables  for  the  five 
variables  in  the  AK  nonlinear  term  for  both  Kw(b)  and  Kb(W)  are  given  in  References  1  and  2  and  will 
not  be  repeated  here. 

As  seen  in  Figure  9,  Kw(B)  in  general  can  deviate  slightly  from  SET  or  LT  near  AOA  0  deg. 
It  then  decreases  until  it  reaches  a  minimum  value  and  then  approaches  a  value  of  1 .0  at  high  AOA. 
On  the  other  hand,  Kb(w)  can  either  increase  or  decrease  past  AOA  0  deg.  Eventually,  it  also 
decreases  until  it  reaches  some  minimum  value  at  high  AOA.  The  physics  of  what  occurs  in  this 
nonlinear  behavior  and  the  details  of  the  interference  factor  nonlinearities  are  given  in  References  1 
and  2. 


The  general  nonlinear  model  for  the  wing-body  interference  due  to  control  deflection  is  of 
the  form 


^W(B)  ^l(^)  [^W(B)lsBT  (29) 

for  both  the  ^  =  0  and  =  45  deg  roll  orientations.  However,  the  values  Ci(M)  and  Cjd  a^l ,  M) 
are  different  for  the  two  roll  positions.  Tables  of  Cj  and  Q  are  given  in  Reference  2.  The  nonlinear 
model  for  the  body-tail  interference  due  to  control  deflection  is  of  the  form 

^B(W)  “  t^B(W)^SBT  (30A) 

for  $  =  0  deg  roll  and  of  the  form 

^B(W)  ~  t^B(W)lsBT  (30E) 

for  ^  =  45  deg  roll.  Again,  the  constants  C(aw)  and  0(6^)  are  defined  in  Reference  2. 

Figure  10  presents  a  qualitative  view  of  how  the  interference  terms  kvv(B)  and  behave. 
In  examining  Figure  lOA,  it  is  seen  that  for  low  Mach  numbers,  slender  body  theory  generally 
underpredicts  the  value  of  kw(B)  at  low  AOA.  At  some  value  of  (%,  which  is  bout  25  deg,  k,y(B)  starts 
decreasing  rapidly  and  goes  to  zero  around  a  value  of  =  50  to  55  deg.  This  behavior  is  believed 
to  be  the  combination  of  wing  stall  and  air  flowing  between  the  wing  and  body  (or  blow-by).  On  the 
other  hand,  for  high  supersonic  Mach  numbers,  just  the  reverse  behavior  occurs.  At  low  AOA,  SET 
overpredicts  the  value  of  kw(B)-  However,  as  AOA  increases,  the  value  of  kw(B)  increases  due  to 
compressibility  nonlinearities  until  a  point  where  the  blow-by  effects  offset  the  compressibility 
effects.  The  addition2il  lift  on  the  wing  then  goes  to  zero  at  around  an  of  75  deg. 
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Figure  lOA  was  derived  primarily  from  total  missile  configuration  data  bases  as  opposed  to 
the  missile  component  data  base  of  Reference  38.  The  reason  for  this  was  the  control  fins  in  the 
Reference  38  data  base  were  too  small  to  get  accurate  values  of  the  additional  lift  due  to  6.  The 
downside  to  the  approach  taken  to  derive  k^(3)  is  the  fact  that  without  missile  force  and  moment 
components  identified,  part  of  the  kw(B)  term  may  actually  be  kB(w).  Hence,  the  model  that  represents 
^Bov)  of  Figure  lOB  basically  represents  the  remaining  load  due  to  control  deflection  required  to 
match  total  missile  configuration  data  bases.  As  seen  in  Figure  lOB,  kB(w)  starts  out  with  its  SBT 
theory  value  at  |  a^l  or  |6w|  =  0  deg  (depending  on  whether  roll  of  0  or  45  deg  is  of  interest),  and 
then  decreases  at  some  value  of  |  or  16^1  in  a  linear  fashion.  CFD  analysis  or  additional  wind 
tunnel  data  with  larger  control  surfaces  would  allow  refinement  of  Equations  (29)  and  (30). 


The  last  interference  term  considered  is  that  of  the  wing-shed  vortices  on  the  tail  surfaces. 
It  is  the  next  to  last  term  of  Equation  (6).  In  general,  the  nonlinear  form  of  Cj^  is  defined  by 


'T(V) 


"N. 


T(V) 


2ti  (AR).j.  -  r^)  Ajgf 


(sj  -  r.j.)  j^ij  cos<^  +  i^sin^j  (31) 


Equation  (31)  reverts  back  to  SBT  for  low  AOA.  At  higher  AOA,  ,  Kw(B),  kw(B),  F,  ij, 

and  all  have  nonlinearities.  The  terms  ij  and  of  Equation  (31)  represent  the  wing-tail 
interference  factor  for  the  windward  and  leeward  planes  respectively.  At  0  =  0  deg,  only  the  ij  factor 
is  of  interest,  whereas  at  ^>  =  45  deg,  both  the  leeward  and  windward  planes  are  considered.  The 
nonlinearity  is  introduced  into  i,  and  i4  in  the  $  =  45  deg  roll  by  defining  loading  factors  for  the 
windward  and  leeward  plane  fins.  These  factors  are: 
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1.0 


1.0  - 


;  a  ^  65 


=  1.6  ,  Pj  =  0.4  ;  a  >  65 


(32) 


Pw  and  P{  have  the  effect  of  increasing  or  decreasing  the  strength  of  the  wing-shed  vortex  from  the 
windward  or  leeward  plane  respectively.  These  factors  are  used  in  the  calculation  of  ij  and  i4. 

The  parameter  F  of  Equation  (31)  is  used  as  a  control  on  the  control  deflection  component 
of  determined  empirically  based  on  comparisons  of  theory  to  experiment.  It  is  defined 

in  the  same  tables  as  kw(B)  and  kB(W)  of  Reference  2.  To  put  an  additional  nonlinear  control  of  the 
first  term  of  Equation  (31)  was  broken  down  into  components  due  to  a  and  6.  That  is. 


(33) 
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The  term  due  to  a,  was  then  defined  as 


Cm  1  =  a  +  Ba  +  Ca^  +  Da^ 


for  ^>  =  0  deg,  and  as 


h  1  =  G,  [Cn  1 


(34) 


(35) 


for  =  45  deg.  A,  B,  C,  D,  and  G,  were  all  constants  determined  empirically  based  on  data  and  are 
also  defined  in  Reference  2.  In  addition,  a  maximum  value  constraint  was  placed  on  the  overall 
value  of  Cm  so  as  not  to  allow  it  to  exceed  the  lift  on  the  tail  alone. 

^T(V) 


3.3  CENTER  OF  PRESSURE 

The  body-alone  center  of  pressure  is  computed  based  on  a  weighted  average  of  the  linear  and 
nonlinear  components  of  normal  force.  That  is 


X. 


^n^(^cp)nl 


CP 


Cn_  +  c 


Nm 


(36) 


where  the  linear  center  of  pressure  is  defined  based  on  the  Section  2  low  AOA  methods  and  the 
nonlinear  center  of  pressure  is  the  centroid  of  the  planform  area.  Two  nonlinear  phenomena  occur 
on  the  body-alone  center  of  pressure  that  Equation  (36)  does  not  fully  account  for.  These  are  the 
asymmetric  shedding  of  vortices  at  AOA  greater  than  about  25  deg  and  at  M„  <  2.0  and  the  transonic 
effects  of  a  shock  standing  on  the  body.  As  a  result,  a  table  of  center  of  pressure  shifts  was  defined 
empirically  as  a  function  of  AOA  and  Mach  number.  This  table  is  given  in  Reference  1. 

The  center  of  pressure  of  the  wing-alone  and  wing-body  normal  force  is  assumed  to  vary 
from  its  linear  value  at  zero  AOA  to  the  centroid  of  the  wing  planform  at  60  deg  AOA.  If  A  and  B 
are  the  centers  of  pressure  of  the  linear  and  nonlinear  normal  force  terms  (in  percent  of  mean 
geometric  chord),  and  =  a  +  6,  then  the  center  of  pressure  of  the  wing-body  or  wing-alone  lift  is 


(XcpU  '  Ww  i“wl  (B-A)  +  ^  4  (A-B)  (37) 


Equation  (37)  is  the  methodology  used  for  roll  position  of  0  deg. 

When  the  fins  are  rolled  to  a  non-zero  roll  orientation,  the  center  of  pressure  Equation  (37) 
will  change  because  of  the  geometry  of  the  wings  and  an  asymmetric  effect  on  the  wing  loading.  To 
visualize  this  effect,  imagine  a  missile  rolled  to  ^>  =  45  deg  and  increasing  in  AOA.  As  AOA 
increases,  two  things  occur.  First,  the  windward  plane  fins  carry  more  and  more  of  the  load 
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compared  to  the  leeward  plane  fins.  Second,  the  local  Mach  number  in  the  windward  plane  is 
different,  and  typically  lower,  than  the  leeward  plane.  This  has  the  effect  of  shifting  the  wing-alone 
center  of  pressure  forward  in  the  windward  plane.  Since  the  load  and  wing  centers  of  pressure  are 
different  on  the  windward  and  leeward  plane  fins,  this  results  in  a  net  forward  shift  in  the  center  of 
pressure  for  O  =  45  deg  roll  compared  to  the  <E>  =  0  deg  computation  of  Equation  (37).  This  shift 
appears  to  occur  for  all  Mach  numbers,  and  is  largest  at  moderate  AOA,  and  goes  to  zero  at  AOA 
0  and  90  deg.  At  90  deg  AOA,  the  windward  plane  fins  carry  almost  all  the  load  compared  to  the 
leeward  plane  fins,  but  geometrically,  the  fins  are  all  aligned  perpendicular  to  the  AOA  plane.  This 
center  of  pressure  shift  was  derived  in  Reference  2  and  is 


r  + 


c. - 


cos($)  sin(2a) 


'  0.8k  ^ 
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;  a  <  65 
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(38) 


Equation  (38)  is  added  to  Equation  (37)  for  the  roll  orientation  of  45  deg. 

The  center  of  pressure  of  the  body-wing  carryover  normal  force  contribution  is  at  the  centroid 
of  the  Mach  box  created  by  the  Mach  lines  from  the  leading  and  trailing  edges  of  the  wings  as  they 
cross  the  body  planar  surface.  For  ^>  =  45  deg  roll,  the  local  Mach  number,  computed  by  Modified 
Newtonian  Theory,  is  used  to  define  the  Mach  box.  This  has  the  effect  of  shifting  the  center  of 
pressure  of  this  component  of  normal  force  slightly  forward. 


3.4  NONLINEAR  AERODYNAMIC  LOADS  DISTRIBUTION 

Section  3.3  discussed  how  the  nonlinear  aerodynamic  loads  were  computed  for  typical 
missile  configurations.  This  nonlinear  aerodynamic  load  is  defined  primarily  in  terms  of  total  force, 
moment  and  center  of  pressure.  Some  distribution  information  is  provided  for  the  nonlinear 
component  of  the  body  loads,  but  none  is  currently  available  for  the  lifting  surface  loads  or  for  their 
interference  effects  on  the  body.  This  situation  is  not  restrictive  in  terms  of  aerodynamic  and 
performance  analysis.  Structural  engineers,  however,  need  to  know  not  only  the  magnitude  of  the 
aerodynamic  forces,  but  also  how  they  are  distributed  over  the  surface  of  a  missile  if  they  are  to 
determine  the  shear  and  bending  moments  to  which  its  components  will  be  subjected.  For  this 
reason,  the  AP95  code  was  modified  to  allow  the  prediction  of  the  distribution  of  the  nonlinear,  as 
well  as  the  linear,  aerodynamic  loads  over  both  the  body  and  control  surfaces.  All  interference 
effects  were  included  in  the  analysis  in  addition  to  the  individual  component  nonlinear 
aerodynamics.  It  is  believed  that  the  next  version  of  the  Aeroprediction  Code  (AP98)  will  be  the  first 
code  of  its  kind  to  offer  this  capability. 
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The  body-alone  load  distribution  is  composed  of  a  linear  term  that  is  computed  based  on 
integration  of  the  pressures  computed  along  and  around  the  body  (see  Section  2)  and  the  local 
component  of  nonlinear  load  from  the  viscous  crossflow  methodology.  That  is. 


/ 


Cp(x,(j))(cos(l))rd4) 


(39) 


Comparison  of  Equation  (39)  results  to  CFD  results^  showed  that  the  predicted  load  in  the  nose 
region  was  slightly  too  high  at  low  AO  A.  As  a  result,  some  of  the  nonlinear  load  in  the  nose  region 
was  shifted  to  the  afterbody  below  AOA  at  30  deg.  The  total  normal  force  and  pitching  moments 
were  held  constant  before  and  after  the  load  redistribution.  A  better  match  of  approximate  body 
loads  with  CFD  was  then  obtained. 

The  body-wing  and  body-tail  loads  of  Equation  (6)  are  represented  using  the  Mach  lines  from 
the  leading  and  trailing  edges  of  the  wing  or  tail  in  conjunction  with  the  total  normal  force 
component  load.  The  total  normal  force  component  load  is  used  to  determine  the  height  of  the 
trapezoidal  load  created  in  the  vicinity  of  the  root  chord  by  the  carryover  of  the  induced  wing  effects 
onto  the  body.  If  these  induced  effects  fall  behind  the  end  of  the  body,  the  height  of  the  trapezoid 
is  readjusted  to  make  sure  the  total  local  integrated  loads  equals  that  of  the  aerodynamic  component. 

The  load  on  the  wing  in  conjunction  with  the  body  of  Equation  (6)  is  treated  in  the  following 
manner.  First,  the  local  linear  load  and  its  center  of  pressure  are  computed  based  on  Section  2 
methods.  The  local  load  along  any  chord  of  the  wing  or  tail  is  then: 


TE 

Ml  '  f  ^C,(x,y)dx  (40) 

^LE 

The  chordwise  center  of  pressure  of  that  load  is  simply  the  center  of  gravity  of  the  local  chordwise 
load.  The  total  load  on  any  wing  chord  is  assumed  to  be  of  the  form 


He  (vll  ) 

C  1 
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r»wjijw 
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(41) 


Equation  (41)  basically  assumes  the  nonlinear  load  is  distributed  similarly  to  the  linear  load.  In 
validating  this  assumption  against  CFD  computations,  it  was  found  to  be  quite  reasonable. 

The  wing-tail  interference  load  (next  to  last  term  of  Equation  (6))  is  assumed  to  be  distributed 
in  a  “l/i^  ”  manner  from  where  the  wing-shed  vortex  hits  the  tail  surface.  Also,  some  of  the  wing-tail 
load  is  distributed  onto  the  body  between  the  tail  surfaces.  The  reader  is  referred  to  Reference  3  for 
the  details  of  this  methodology. 
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When  the  missile  is  in  the  roll  position  of  O  =  45  deg,  the  body  carryover  load  shows 
somewhat  of  a  different  behavior  than  at  <E>  =  0  deg.  At  =  45  deg,  the  high  pressure  region  in  the 
vicinity  of  the  wings  is  maintained  just  like  at  =  0  deg.  However,  behind  the  wings,  the  body 
carryover  can  actually  go  negative.  This  is  modeled  by  allowing  the  positive  trapezoidal  load 
induced  by  the  Mach  lines  from  the  wing  leading  and  trailing  edges  to  go  negative  behind  the  root 
chord.  The  only  constraint  is  that  the  total  load  and  pitching  moment  remain  constant  at  the  values 
determined  by  the  aerodynamic  load  and  pitching  moments.  Also,  the  windward  plane  wing  carries 
more  load  than  the  leeward  plane  wing.  This  is  modeled  through  the  use  of  Equation  (32). 


4.0  ASYMMETRIC  CONFIGURATION  METHODOLOGY 


The  body-alone  nonaxisymmetric  body  methodology  was  based  on  extensions  of  the  method 
of  Jorgensen,^^  and  the  wing-body  interference  effects  methodology  was  based  on  extensions  of 
Nelson"*^  and  Est  and  Nelson^’  at  low  AOA  and  the  method  of  the  present  authors^’^  at  high  AOA. 
Several  new  additions  to  the  state-of-the-art  were  made  in  extending  the  methods  of  References  45 
through  47.  Figure  2  shows  the  noncircular  wing-body  cases  of  interest. 

Basically,  the  method  of  Jorgensen,  extended  by  the  Reference  5  methods,  defined  the  axial 
and  normal  force  and  pitching  moment  coefficients  as 


Ca  =  Ca^  cos^a 


(42) 
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The  subscripts  L  and  NL  represent  the  linear  and  nonlinear  components  of  normal  force  and  pitching 
moment  respectively.  The  factors  (C„/C„  )sb  and  represent  the  slender-body  and 

Newtonian  approximations  to  the  local  normtS  force  coefficient  per  unit  length  of  the  desired  cross- 
sectional  shape  (C„)  to  the  similar  coefficient  for  the  equivalent  circular  cross-sectional  shape  (C  ) . 
The  major  contribution  of  Jorgensen  was  in  showing  that  aerodynamics  of  a  noncircular  cross- 
sectional  shape  could  be  computed  reasonably  accurately  to  high  AOA  if  the  factors  C  /C  were 
defined.  He  defined  these  for  ellipses.  " 
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Jorgensen’s  theory  was  more  limited  than  desired  to  consider  all  the  cases  of  Figure  2.  As 
a  result.  Reference  5  provided  significant  extensions  to  the  Jorgensen  methodology.  First  of  all, 
additional  equations  were  derived  for  (Cj^/C^j  )j^  for  shapes  other  than  ellipses  (see  Figure  1 1). 
Secondly,  Newtonian  correction  factors  (NF)  were  derived  for  low  crossflow  Mach  numbers  where 
Newtonian  theory  was  inaccurate.  Jorgensen  had  no  factor,  NF,  in  his  equations.  Figures  12  through 
14  give  these  values  along  with  the  equations  used  to  compute  NF  for  ellipses,  squares,  and  triangles. 
Third,  a  critical  Reynolds  number  was  defined  for  noncircular  shapes  as  a  function  of  the  circular 
value  of  Rj^  .  These  values  are  given  in  Figure  15.  Fourth,  the  AO  A  axial  force  was  computed 
using  References  1  and  4  for  the  equivalent  axisymmetric  body.  The  only  other  change  in  axial  force 
coefficient  was  in  the  skin-friction  component  where 


(^A,lc  ■  Nc 


(Cir)Nc 

(Cir)c 


(45) 


Fifth,  wing-body  interference  factors  were  derived  based  on  the  methods  of  References  46  and  47 
at  low  AOA  (see  Figure  16).  At  high  AO  A,  wing-body  interference  factors  were  derived  based  on 
the  methods  of  References  1  and  2.  Thus,  the  body  cross-sectional  shape  was  allowed  to  influence 
the  aerodynamics  at  low  and  moderate  AOA,  whereas  at  AOA  greater  than  about  30  deg,  it  was  not. 
Finally,  slender  body  theory  scaling  factors  were  derived  to  allow  one  to  compare  aerodynamics 
computed  based  on  an  equivalent  axis5nnmetric  body  or  a  body  of  constant  width  to  wind  tunnel  data 
taken  on  a  body  of  constant  width  or  constant  cross-sectional  area.  Figure  17  gives  a  summary  of 
the  slender  body  theory  scaling  factors  for  use  in  Figure  16B.  Figure  16A  was  based  on  7tdgqV4  and 
all  elliptical  wind  tunnel  data  examined  by  the  author  was  taken  that  way,  so  (SBTSF),  and  (SBTSFlj 
of  Figure  17  are  both  1. 

The  brief  summary  given  here  does  not  go  into  any  of  the  details  of  the  derivation  of 
Figures  12  through  17  or  the  equations  that  accompany  them.  The  interested  reader  is  once  again 
referred  to  Reference  5  for  these  details. 

One  final  point  needs  to  be  made  before  leaving  the  summary  on  nonaxisymmetric  body 
aerodynamics.  That  concerns  the  local  loads  used  by  the  structural  analyst.  Reference  3  developed 
a  fairly  accurate  way  of  distributing  the  nonlinear  aerodynamic  loads  over  the  body  and  lifting 
surfaces  based  on  CFD  results  for  an  axisymmetric  body.  What  is  generated  within  the 
aeroprediction  code  is  the  aerodjmamic  load  distribution  for  the  equivalent  circular  body.  It  is 
possible  that  some  noncircular  cross  sections  may  have  loading  patterns  that  are  significantly 
different.  Unfortunately,  there  was  not  enough  data  available  to  determine  whether  this  was  the  case 
or  to  attempt  to  model  these  effects  if  present.  To  resolve  this  issue  would  require  either  wind  tunnel 
measurements  or  CFD  computations  to  determine  detailed  surface  pressure  maps  for  an  extensive 
range  of  cross  section  shapes.  Such  an  endeavor  was  beyond  the  scope  of  the  present  work,  but 
should  certainly  be  considered  in  the  future  if  sufficient  time  and  resources  became  available. 
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b.  Squares 


HGURE  11.  RATIO  OF  LOCAL  NORMAL  FORCE  COEFHCIENT  FOR  VARIOUS  CROSS  SECTIONAL 
BODY  SHAPES  TO  THAT  FOR  THE  EQUIVALENT  CIRCULAR  CROSS  SECTION 
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Low  AOA  Merging 
NF  =  1  +  [(NF),  -1]  —  ;  a  ^  20 


NF  =  (NF)j  ;  a  >  20 


Values  of  a/b  other  than  0.5.  2.0.  3.0 


NF= 

;  a/b>3.0 

NF=  {(NF)^=3  -  (NFU=2}  (a^  -  2)  .  (m^.^ 

;  2^a/b^3.0 

NF=  {(NF)3;b=2  -  1}  (a^  -1)  +  1 

;  l^a/b<2.0 

NF=  {(NF)^=o.5  -  1}  2(1  -  a^)  +  1 

;  0.5^a/b<1.0 

(l'lP)a/b=0.5 

;  a/b^O.5 

HGURE  12.  NEWTONIAN  CORRECTION  FACTOR  METHODOLOGY  FOR  AN  ELLIPSE 
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Newtonian  Correction  Factor  for  Squares  as 
a  Function  of  Comer  Radius  at  =  0 


Newtonian  Correction  Factor  for  Squares  as 
a  Function  of  Crossflow  Mach  Number 
(k  =  0) 
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[(NF)m^=o.55  -  1]  ;  0.55  <  ^  2.6 


>  2.6 


=  0.4  -  2.67k  ;  k  s  0.15 
=  0  ;  k  >  0.15 

FIGURE  13.  NEWTONIAN  CORRECTION  FACTOR  METHODOLOGY  FOR  SQUARES 
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Newtonian  Correction  Factor  for  Triangles 
as  a  Function  of  Comer  Radius  (M^  =  0) 


Triangle 
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Mn  - 

NF=  NF„  - 
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Newtonian  Correction  Factor  for  Triangles 
as  a  Function  of  Crossflow  Mach  Number 
(k  =  0) 
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Inverted  Triangle 
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1.0  ;  Mn>2.4 


M„  =  0.4  -  2.67k  ;  k  ^  0.15 
M^,  =  0  ;  k  >  0.15 

HGURE  14.  NEWTONIAN  CORRECTION  FACTOR  METHODOLOGY  FOR  TRIANGLES 
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k 

Critical  Reynolds  Number  for  Squares  and  Triangles 


0  12  3 

a/b 

Critical  Reynolds  Number  for  Ellipses 


where  Rj^^  Defined  for  a  given  W/T  or  Flight  Test 
R„  Typically  1  x  10*  to  4  x  10* 

FIGURE  15.  CRITICAL  REYNOLDS  NUMBER  FOR  NONCIRCULAR  SHAPES 
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Impact  of  Square  and  Triangle  Body  Cross  Sections  on  Low  Mach  Number  Values  of  K^(B) 
(Data  from  References  51  and  52  And  Based  on  =  n  WV4) 


HGURE  16.  WING-BODY  INTERFERENCE  FACTORS  (LOW  AOA) 
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Note:  Equations  (1)  through  (4)  apply  to  triangles  and  squares  of  Figure  16B.  Ellipses  (Figure  16A)  derived 

based  on  izd^^/4.  (s/r),  is  local  value  of  s/r  where  wing  is  located. 


HGURE  17.  SBT  SCALING 
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5.0  SUMMARY  OF  AERODYNAMIC  METHODS 


A  summary  of  the  theoretical  and  empirical  methods  used  for  computing  aerodynamics  in 
the  AP98  are  given  in  Figures  18  through  20.  Figure  18  summarizes  the  body-alone  methods, 
Figure  19  summarizes  the  wing  and  interference  aerodynamic  methods,  and  Figure  20  summarizes 
the  dynamic  derivative  methods.  All  the  methods  shown  in  Figures  18  through  20  were  discussed 
in  Sections  2  through  4  of  this  report.  The  major  new  additions  to  the  AP98  are  shown  in  Figures  18 
and  19,  which  include  the  Improved  Axial  Force  at  AO  A,  ^>  =  45  deg  aerodynamics.  Nonlinear 
Structural  Loads  and  Nonaxisymmetric  Body  Aerodynamics.  No  changes  have  been  made  to  the 
d5mamic  derivative  methods  of  Figure  20,  which  are  still  linear  with  no  nonlinear  terms  included. 
Also,  although  not  shown  in  Figures  18  through  20,  a  significantly  more  robust  pre-  and  post¬ 
processing  software  package  for  personal  computer  usage  will  be  developed  for  the  AP98.  Using 
a  personal  computer  with  a  200  megahertz  chip,  a  single  case  can  be  executed  in  much  less  than  a 
second,  even  though  the  computer  code  listing  now  exceeds  17,000  lines. 


COMPONENT/ MACH 
NUMBER  REGION 

SUBSONIC 

M,  <  0.8 

TRANSONIC 

0.8  ^  M^  :s  1.2 

LOW 

SUPERSONIC 

1.2  ^  M^  1.8 

MOD/HIGH 

SUPERSONIC 

1.8  ^  M,  6.0 

HYPERSONIC 

M,  >  6.0 

NOSE  WAVE  DRAG 

EMPIRICAL 
(Ref.  6) 

SEMEMPIRICAL 
BASED  ON  EULER 
SOLUTIONS 
(Ref.  19) 

SECOND-ORDER 
VAN  DYKE  PLUS 
MNT  (Ref.  6) 

SOSETPLUS  IMNT 
(Ref.  13) 

SOSETPLUS  IMNT 
MODEED  FOR 
REAL  GASES 
(Ref.  14) 

BOATTAE.  OR  FLARE 
WAVE  DRAG 

— 

WUAND 
AOYOMA 
(Ref.  6) 

SECOND-ORDER 
VANDYKE 
(Ref.  7) 

SOSET(Ref.  11) 

SOSET  FOR  REAL 
GASES  (Ref.  14) 

SKIN  FRICTION  DRAG 

VAN  DRIEST  E 
(Ref. 15) 

BASE  DRAG 

IMPROVED  EMPIRICAL  METHOD 
(Ref. 17) 

AXIAL  FORCE  AT  a 

IMPROVED  EMPIRICAL  METHOD 
(Ref.  4) 

AEROHEATING 

INFORMATION 

— 

SOSET  PLUS  IMNT  FOR  REAL  GASES 
(Ref.  21) 

INVISCID  LIFT  AND 
PITCHING  MOMENT 

EMPIRICAL 
(Ref.  6) 

SEMEMPIRICAL 
BASED  ON  EULER 
SOLUTIONS 
(Ref.  19) 

TSEN  FIRST- 
ORDER 
CROSSFLOW 
(Ref.  8) 

SOSET  (Ref.  11) 

SOSET  FOR  REAL 
GASES  (Ref.  14) 

VISCOUS  LIFT  AND 
PITCHING  MOMENT 

IMPROVED  ALLEN  AND  PERKINS  CROSSFLOW 
(Ref.  42) 

NONAXISYMMETRIC 
BODY  AERO 
(^>  =  0,45°) 

MODIFIED  JORGENSEN 
(Ref.  5) 

NONLINEAR  ST. 
LOADS  AVAIL. 

(^>  =  0,45°) 

NO 

YES 
(Ref.  3) 

HGURE  18.  AP98  METHODS  FOR  BODY-ALONE  AERODYNAMICS 


43 


NSWCDD/TR-98/1 


COMPONENT/ MACH 
NUMBER  REGION 

SUBSONIC 

M^  <  0.8 

TRANSONIC 

0.8  £  M^  s  1.2 

LOW 

SUPERSONIC 

1.2  ^  M,  ^  1.8 

MOD/HIGH 

SUPERSONIC 

1.8  ^  ^  6.0 

HYPERSONIC 
M^  >  6.0 

WAVE  DRAG 

EMPIRICAL 
(Ref.  27) 

LINEAR  THEORY 
PLUS  MNT 
(Ref.  27) 

SHOCK 

EXPANSION  (SE) 
PLUS  MNT  ALONG 
STRIPS  (Ref.  13) 

SE  PLUS  MNT 
FOR  REAL 
GASES  ALONG 
STRIPS  (Ref.  13) 

SKIN  FRICTION  DRAG 

1  VAN  DRIEST  n  (Ref.  15) 

TRAILING  EDGE 
SEPARATION  DRAG 

EMPIRICAL  (Ref.  27) 

BODY  BASE  PRESSURE 
CAUSED  BY  TAtt,  FINS 

IMPROVED  EMPIRICAL  (Ref.  17) 

INVISCID  LIFT  AND 
PITCHING  MOMENT 

•  LINEAR 

•  LIFTING 
SURFACE 
THEORY  (Ref.  27) 

•  EMPIRICAL 
(Ref.  27) 

•  3DTWT 
(Ref.  27) 

•  3DTWTORSE 
(Ref.  27  or  13) 

•  3DTWTORSE 
(Ref  27  or  14) 

•  NONLINEAR 

•  EMPIRICAL  (Ref.  1,  2) 

WING-BODY,  BODY¬ 
WING  INTERFERENCE 
(^>  =  0, 45*’) 

•LINEAR 

•  SLENDER  BODY  THEORY  OR  LINEAR  THEORY  MODIFIED 

FOR  SHORT  AFTERBODIES  (Ref.  1, 2) 

•NONLINEAR 

•  EMPIRICAL  (Ref.  1,2) 

WING-BODY, 
INTERFERENCE  DUE  TO  5 
($  =  0,45°) 

•  LINEAR 

•  SLENDER  BODY  THEORY  (Ref  2) 

•  NONLINEAR 

•  EMPIRICAL  (Ref.  1,2) 

WING-TAIL 
INTERFERENCE 
(^  =  0,45") 

UNE  VORTEX  THEORY  WITH  MODIHCATIONS  FOR  Kw®, 

TERM  AND  NONUNEARITIES  (Ref  2) 

AEROHEATING 

NONE  PRESENT 

SE  PLUS  MNT 
(Ref.  21) 

SE  PLUS  MNT 
REAL  GASES 
(Ref.  21) 

NONAXISYMMETRIC 
BODY  AERO  ($  =  0,45°) 

IMPROVED  NELSON  ESTIMATE  FOR  AP98  (Ref  5, 46, 47) 

NONLINEAR  ST.  LOADS 
AVAIL.  ($  =  0,45°) 

NO 

YES  (Ref.  3) 

FIGURE  19.  AP98  METHODS  FOR  WING-ALONE  AND  INTERFERENCE  AERODYNAMICS 
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COMPONENT/ 
MACH  NUMBER 
REGION 

SUBSONIC 

M.  <  0.8 

TRANSONIC 

0.8  ^  M^  <  1.2 

LOW 

SUPERSONIC 

1.2  i  M,  5  1,8 

MOD/HIGH 

SUPERSONIC 

1.8  ^  ^  6.0 

HYPERSONIC 
>  6.0 

BODY  ALONE 

EMPIRICAL 

WING  AND 
INTERFERENCE 
ROLL  DAMPING 
MOMENT 

LIFTING 

SURFACE 

THEORY 

EMPIRICAL 

LINEAR  THIN 
WING  THEORY 

LINEAR  THIN  WING  OR 

STRIP  THEORY 

WING  MAGNUS 
MOMENT 

ASSUMED  ZERO 

WING  AND 
INTERFERENCE 
PITCH  DAMPING 
MOMENT 

LIFTING 

SURFACE 

THEORY 

EMPIRICAL 

LINEAR  THIN 
WING  THEORY 

LINEAR  THIN  WING  OR 

STRIP  THEORY 

HGURE  20.  AP98  METHODS  FOR  DYNAMIC  DERIVATIVES  (REFERENCES  28  AND  16) 


6.0  RESULTS  AND  DISCUSSION 


Several  cases  are  considered  to  try  to  compare  the  new  AP98  with  experimental  data,  the 
AP95,  and  other  theoretical  methods  (if  computations  are  available).  In  general,  it  is  believed  the 
AP98  comparisons  to  experimental  data  are  as  good  as  or  better  than  the  AP95.  However,  the  AP98 
is  more  robust  in  terms  of  roll  orientation  and  configuration  cross  sectional  geometry.  Hence,  for 
<})  =  45  deg  aerodynamics,  only  AP98  results  will  be  shown.  Also,  for  nonaxisymmetric  bodies,  only 
AP98  results  will  be  shown. 


6. 1  AXIS  YMMETRIC  BODY  CONHOURATION  CASES 

The  first  axisymmetric  body  case  considered  is  a  canard-body-tail  configuration  shown  in 
Figure  21  A.  This  case  was  tested^^  at  M„  =  0.2  at  the  Naval  Postgraduate  School  in  Monterey.  The 
configuration  tested  had  aspect  ratio  canards  of  1.59  and  tails  of  0.9.  The  body  was  22.6  calibers  in 
length.  Test  data  were  available  at  50  deg  AOA  for  both  4>  =  0  and  45  deg  roll  orientation  with 
canards  deflected  0  and  ±20  deg.  Reference  53  pointed  out  that  the  sting  balance  was  designed  for 
the  large  normal  forces  at  high  AOA  and  hence  the  axial  force  measurements  are  not  expected  to  be 
as  accurate.  The  wind  tunnel  model  with  boundary  layer  trip  option  was  chosen  for  the 
aeroprediction  computations.  Sea  level  conditions  were  also  assumed.  Figures  21B,  21C,  and  21D 
give  the  axial  and  normal  force  and  pitching  moment  coefficients  for  the  roll  position  of  $  =  0  and 
control  deflections  of  0  deg  and  ±20  deg.  Results  are  given  for  the  experimental  data,  AP95,  AP98 
and  Missile  DATCOM^  (which  was  taken  from  Reference  53).  As  seen  in  the  figures,  the  AP95  and 
AP98  give  similar  results  for  normal  force  and  pitching  moments  and  are  both  in  reasonable 
agreement  with  data.  The  AP98  gives  some  slight  improvement  over  the  AP95  in  axial  force  at 
higher  AOA  due  to  the  new  axial  force  methodology  included.  Both  the  AP95  and  AP98  appear  to 
give  better  results  compared  to  data  than  the  DATCOM  results  presented  in  Reference  53. 
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CANARD 


NOTE:  ALL  DIMENSIONS  IN  INCHES,  FULL  SCALE 
FULL 

TAIL  \  ^  K-  2.0 


98.3“ 


113.0 


t 

5.0 


NOSE  RADIUS  1.42 


>1  |<-0.75  ->|  1<-1.06 


HGURE  21  A.  CANARD-CONTROLLED  MISSILE  CONHGURATION"^ 
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FIGURE  21 C.  COMPARISON  OF  STATIC  AERODYNAMICS  BETWEEN  EXPERIMENT  AND 
THEORY  FOR  HGURE  21 A  CONHGURATION  (6  =  0  DEG,  $  =  0  DEG,  M„  =  0.2) 
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The  second  axisymmetric  body  case  is  taken  from  Reference  55.  This  data  set  consisted  of 
various  length  to  diameter  canard-body-tail  configurations  tested  at  =  2.01  to  AO  A  25  to  30  deg. 
Data  were  taken  at  a  R^/ft  of  3.47  x  10®  and  with  no  boundary  layer  trip  present.  Figure  22 A  shows 
the  configurations  considered  for  comparison  to  data.  These  cases  consist  of  19.1,  16.7  and 
14.8  caliber  canard-body-tail  configurations.  The  tail  has  about  2  percent  of  the  tip  eliminated,  but 
this  was  not  considered  in  the  AP95  and  AP98  results.  Hence,  it  should  be  expected  that  the  theory 
should  be  slightly  high  on  normal  force  and  slightly  more  stable  pitching  moments  compared  to  data 
Figure  22B  shows  the  normal  force  and  pitching  moment  coefficient  comparisons  of  the  AP95  and 
AP98  to  data.  Even  with  the  slightly  larger  tail  fins  used  in  the  theory,  the  normal  force  and  center 
of  pressure  average  errors  are  well  within  the  ±10  percent  and  ±4  percent  of  body  length  goals 
respectively.  The  AP98  normal  force  predictions  are  slightly  better  than  the  AP95,  whereas  the 
opposite  is  true  of  the  pitching  moments. 


Vertical  5.4S0 


HGURE  22A.  CANARD-BODY-TAIL  CONHGURATIONS  OF  VARIOUS  LENGTH  TO  DIAMETER 
RATIOS  (REFERENCE  55)  (ALL  DIMENSIONS  IN  INCHES) 
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The  third  case  considered  is  taken  from  Reference  56  and  is  shown  in  Figure  23.  The  model 
was  about  22  calibers  in  length  with  a  sharp  nose  of  2.25  calibers.  The  canards  had  an  aspect  and 
taper  ratio  of  2.0  and  0.3  respectively.  Various  tail  fin  spans  were  considered.  This  model  was 
tested  at  Mach  numbers  1.6  to  3.5  at  AOA  to  about  18  to  20  deg.  It  had  a  boundary  layer  trip  present 
and  was  tested  at  a  R^/ft  of  2.0  x  10^.  Reference  56  gave  separate  values  of  base  axial  force 
coefficient,  which  were  added  to  the  axial  force  values  given  in  the  reference  to  compare  to  the  AP98 
computations.  To  compare  the  experimental  data  to  theory,  Mach  numbers  of  2.5  and  3.5  are 
selected  at  roll  angle  45  deg.  Also,  values  of  the  tail-to-canard  semispan  of  0.47  and  1.25  are 
considered.  Since  0  =  45  deg  roll  is  chosen,  no  AP95  computations  will  be  shown.  Figure  24 
presents  the  comparison  of  theory  to  experiment  for  b/b^.  =  0.47  and  b/b^  =  1 .25.  Results  are  shown 
for  C^,  Cn  and  C^.  As  seen  in  Figure  24,  comparison  of  theory  to  experiment  is  quite  acceptable  and 
meets  the  average  accuracy  goal  of  ±10  percent  on  axial  and  normal  force  and  ±4  percent  of  body 
length  for  center  of  pressure.  The  worst  case  error  on  pitching  moment  at  M  =  2.5  and  b/b^  =  1 .25 
represents  an  error  of  less  than  3  percent  of  the  body  length  in  terms  of  center  of  pressure.  Also 
shown  on  Figure  24  are  M3HAX  computations  taken  from  Reference  57  for  the  b/b^  =  1.25  case  at 
M  =  2.5.  The  results  for  =  0  deg  from  Reference  57  were  assumed  to  apply  to  the  $  =  45  deg 
case  since  no  control  deflections  were  assumed  for  this  example.  As  seen  in  the  figure,  M3HAX 
gives  reasonable  comparison  to  data  at  O  =  45  deg.  Reference  57  shows  M3HAX  giving  improved 
comparisons  to  data  at  $  =  0  deg  roll,  compared  to  ^  =  45  deg  roll,  particularly  for  normal  force  and 
pitching  moment  coefficients. 


The  next  case  considered  is  a  wing-body-tail  case  with  a  low  aspect  ratio  wing  and  tested  at 
the  Naval  Postgraduate  SchooF*  at  M„  =  0.1.  The  configuration  is  shown  at  the  top  of  Figure  25. 
Note  that  the  model  tested  in  the  wind  tunnel  at  the  top  of  the  figure  is  slightly  different  than  that 
where  calculations  were  performed  with  the  AP95  and  AP98.  This  is  due  to  the  fact  the  APC  cannot 
handle  the  detailed  dorsal  and  wing  geometry  shown  at  the  top  of  the  figure.  The  configuration 
shown  in  the  middle  of  Figure  25  has  the  same  aspect  ratio,  span,  taper  ratio,  leading  edge  sweep 
angle  and  area  of  both  the  dorsal  and  tail  as  the  actual  model  at  the  top  of  Figure  25.  Also,  although 
not  shown,  the  distance  to  the  centroid  of  the  planform  areas  is  also  held  constant.  Of  course,  the 
body  is  also  the  same  between  the  two  cases. 

The  only  results  given  in  Reference  58  were  for  normal  force.  Experimental  data  and  Missile 
DATCOM^  results  were  both  given  in  Reference  58  for  roll  positions  of  =  0  and  45  deg.  These 
results  are  compared  to  the  AP95  and  AP98  at  the  bottom  of  Figure  25  for  <I»  =  0  deg  and  to  the 
AI^8  only  for  =  45  deg.  Note  that  comparisons  to  data  are  reasonable  and  quite  good  except  at 
AOA  40  to  70  deg.  It  is  possible  that  wind  tunnel  sting  to  model  interference  increased  the 
experimental  results  in  this  range  somewhat.  Previous  wind  tunnel  studies^®’^  have  concluded  the 
model  sting  can  increase  model  normal  force  loads  by  10  to  20  percent  in  the  high  AOA  range.  Both 
the  AP95  and  AP98  give  slightly  better  comparisons  to  data  than  the  Missile  DATCOM®'‘  for  this 
case.  It  is  suspected  the  low  aspect  ratio  lifting  surfaces,  which  have  a  great  deal  of  nonlinear  lift, 
is  the  primary  reason  for  the  Missile  DATCOM  accuracy  problems,  particularly  at  #  =  0  deg. 
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HGURE  23.  CANARD-BODY-TAIL  CONnOURAHON  WITH  VARYING  TAIL  SPAN 

(ALL  DIMENSIONS  IN  INCHES) 
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HGURE  25.  NORMAL  FORCE  COEFFICIENT  COMPARISONS  OF  BODY-DORSAL-TAIL  (M„  =  0.1) 
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The  next  axis)anmetric-body  missile  configuration  chosen  is  a  wing-body-tail  case  shown 
in  Figure  26,  where  the  wings  are  used  for  control.  Experimental  results  for  this  case  were  taken 
from  Reference  61.  This  configuration  has  a  length  of  about  18  calibers  with  a  tangent  ogive  nose 
2.25  calibers  in  length.  It  has  wings  and  tails  of  fairly  high  aspect  ratios  of  2.8  and  2.6  respectively. 
Data  were  taken  at  Mach  numbers  1.5  to  4.63,  for  AO  As  to  45  deg  and  control  deflections  of  0  and 
10  deg  at  M  of  1 .5  and  2.0  and  0  to  20  deg  at  M  of  2.35  to  4.63.  The  data  were  taken  at  a  Reynolds 
number  of  2.5  x  lOVft  and  boundary  layer  trips  were  also  used.  The  model  had  a  hollow  chamber, 
and  chamber  axial  force  measurements  were  given  separately  in  Reference  20.  These  results  were 
added  to  the  forebody  axial  force  measurements  to  compare  with  the  AP95  and  new  AOA  axial  force 
prediction  method  presented  in  this  report. 

Figure  27  shows  the  comparisons  of  the  AP98  and  AP95  to  the  data  of  Reference  61  for 
O  =  0  deg  and  $  =  45  deg.  Figure  27A  and  27B  give  Q,  Cn  and  for  M  =  1.5  at  6  =  0  and 
6w  =  10  deg  at  ^>  =  0  deg.  In  general,  both  the  AP95  and  AP98  give  acceptable  comparisons  to  data, 
with  the  AP98  giving  slightly  better  results  on  axial  force  at  AOA  than  the  AP95.  On  the  other  hand, 
the  pitching  moments  of  the  AP95  are  slightly  better  than  those  of  the  AP98.  Figure  27C  and  27D 
give  similar  results  for  M  =  2.87,  and  Figure  27E  and  27F,  for  M  =  4.6.  Overall,  for  this 
configuration,  at  =  0  deg  roll,  the  AP95  and  AP98  are  about  equal  in  overall  accuracy 
comparisons.  The  worst  case  errors  are  for  center  of  pressure  at  higher  Mach  number  and  AOA, 
where  the  bow  shock  intersects  the  wing  shocks.  This  nonlinear  phenomena  is  not  modeled  in  the 
=  0  deg  roll  orientation  at  all.  For  the  =  45  deg  roll,  the  center  of  pressure  shift.  Equation  (38), 
partially  accounts  for  this  phenomena,  but  not  entirely.  Center  of  pressure  errors  approach  a  full 
caliber  or  5  percent  of  the  body  length  at  M  =  4.6  and  a  =  40  deg. 

Figure  27G  and  27H  present  M  =  1.5  results  for  =  45  deg  roll.  Since  the  AP95  is  not 
applicable  to  45  deg  roll,  only  the  AP98  and  experimental  data  are  available.  Figure  27G  and  27H 
give  Cn  and  C^  for  M  =  1.5  and  6  =  0  and  10  deg.  Figure  271  and  27J  give  similar  results  for 
M  =  2.87  and  6  =  0  and  20  deg,  and  Figure  27K  and  27L  give  results  for  M  =  4.6  and  6  =  0  and 
20  deg.  Note  that  good  agreement  is  obtained  between  experimental  data  and  the  AP98  for  all  static 
aerodynamics  at  all  three  Mach  numbers  and  for  all  control  deflections.  Here,  the  worst  case  center 
of  pressure  error  is  less  than  3  percent  of  the  body  length. 

Figure  28  presents  results  for  the  Figure  26  configuration  where  the  tail,  versus  the  wing  is 
used  for  control.  Here,  only  the  6^  =  - 10  deg  and  -20  deg  results  are  shown  as  the  6^  =  0  deg 
control  is  basically  the  same  as  that  in  Figure  27  for  6w  =  0  deg.  Both  the  AP95  and  AP98  give 
acceptable  results  for  the  $  =  0  deg  roll  for  C^  and  Cj^^.  However,  the  AP98  is  clearly  superior  to 
the  AP95  for  axial  force  at  high  AOA  and  Mach  number.  Figures  28B,  D  and  F  give  the  tail  control, 
$  =  45  deg  roll  results.  Again,  only  the  AP98  and  experimental  results  are  given.  Note  the  good 
agreement  of  the  AP98  with  the  data.  Again,  the  worst  case  error  on  pitching  moment  results  is  a 
center  of  pressure  error  of  less  than  3  percent  of  the  body  length. 

In  viewing  the  comparison  of  theory  to  experiment  in  Figures  27  and  28  for  both  the  wing 
and  tail  control  alternatives  of  the  Figure  26  configuration,  several  conclusions  can  be  drawn.  In 
terms  of  accuracy,  both  the  AP95  and  AP98  give  acceptable  predictions  for  C^,  and  Cm  for 
=  0  deg  roll  when  a  and  6  are  of  the  same  sign.  The  AP98  gives  superior  results  to  the  AP95  for 
Ca  when  a  and  6  are  of  opposite  signs.  The  AP98  gives  equally  good  results  for  ^  =  45  deg  roll  as 
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for  $  =  0  deg  roll,  whereas  the  AP95  is  only  applicable  for  O  =  0  deg  roll.  The  AP98  gives  slightly 
better  center  of  pressure  predictions  than  the  AP95,  primarily  due  to  Equation  38. 


HGURE  26.  WING-BODY-TAIL  CONHGURATION  USED  IN  VALIDATION  PROCESS** 
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The  next  case  considered  in  the  validation  process  is  shown  in  Figure  29.  This  configuration 
was  tested®^  with  a  body  alone,  body  with  flare  and  with  fins.  The  options  chosen  here  for 
comparison  are  the  body  alone  and  body  with  a  10  deg  flare  and  body  with  15  deg  fins.  The  basic 
body  of  Figure  29  consists  of  a  5-caliber  ogive  forebody  with  a  20  percent  blunt  nose.  The  fins  were 
flat  plates  having  rounding  leading  edges  and  blunt  trailing  edges.  Tests  were  conducted  at  M  =  2.01 
and  Rj/ft  of  2  x  10®.  Base  axial  force  measurements  were  subtracted  out  of  the  experimental  data 
so  the  axial  force  will  not  include  the  base  drag  term.  No  boundary  layer  trip  was  mentioned  in 
Reference  62  so  the  “smooth  model  with  no  boundary  layer  trip”  option  was  used  in  both  the  AP95 
and  AP98  computations.  Data  were  given  to  AOA  24  deg  to  28  deg. 

Figure  30  presents  the  AP95  and  AP98  computations  of  ,  Cj,  and  compared 

to  the  data  of  Reference  62.  Comparison  of  the  theory  to  data  is  acceptable.  However,  the 
comparisons  are  not  as  good  as  desired.  It  is  suspected  that  the  process  of  subtracting 
out  the  base  pressure  from  experiment  and  theory  is  the  primary  source  of  the  disagreement.  To  get 
accurate  base  pressure  measurements  from  experimental  data  generally  requires  at  least  a  three  to 
one  ratio  of  base  to  sting  diameter  plus  several  base  pressure  taps.  Also  the  values  of  wave  plus  skin 
friction  drag  given  by  the  experiment  for  the  body  alone  case  of  0.065  to  0.085  appear  low.  If  the 
base  pressure  component  were  too  high  by  0.02  to  0.04,  this  could  account  for  the  error  between 
experiment  and  theory  for  the  axial  force  coefficient  on  Figure  30A.  Figure  30B  gives  the  same 
static  aerodynamic  results  for  the  10  deg  flare  case  of  Figure  29;  and  Figure  30C,  the  results  for  the 
15  deg  fin  case  of  Figure  29.  Comparison  of  the  AP98  and  AP95  to  experiment  is  excellent,  with 
the  exception  of  the  base  pressure  issue  discussed  previously. 

The  next  set  of  wind  tunnel  data  considered  for  comparison  purposes  is  taken  from 
Reference  63.  Body  alone,  body-tail  and  wing-body-tail  configurations  were  all  a  part  of  this  test 
series.  Figure  31 A  shows  one  of  the  configurations  tested  and  considered  here  for  validation  of  the 
AP98  results.  The  model  is  13.5  calibers  in  length  with  a  1.5  caliber  tangent  ogive  nose.  The  wing 
surfaces  are  fairly  large  with  thickness  of  t/c,  =  0.0178  and  wedge  angles  on  the  leading  and  trailing 
edges  of  15  deg.  The  tail  surfaces  have  thickness  of  t/c^  =  0.05  and  wedge  angles  of  20  deg.  The 
tests  were  conducted  at  Mach  numbers  0.7  to  3.08  with  Reynolds  number  varying  from  about  2  x  10® 
to  4.6  X  10®  per  foot.  The  smooth  model  without  boundary  layer  trip  option  was  used  for  the  AP98 
calculations.  AOA  to  25  deg  were  considered  in  the  wind  tunnel  test.  For  comparison  purposes, 
normal  force  and  pitching  moments  are  compared  to  data  at  M  =  1.42  and  3.08  for  the  =  0  deg  roll 
orientation.  Figure  3  IB  presents  these  results.  As  seen  in  the  figure,  both  the  AP95  and  AP98  give 
quite  acceptable  comparisons  to  data,  with  the  AP98  giving  slightly  better  comparisons  for  pitching 
moment.  The  AP95  gives  slightly  better  normal  force  coefficient  comparisons  to  data  at  the 
M  =  3.08  case  for  a  >  18  deg.  Both  versions  of  the  APC  give  aerodynamics  well  within  the  accuracy 
goals.  Reference  63  also  gave  axial  force  information  where  the  base  pressure  had  been  subtracted 
out.  Unfortunately,  only  a  side  chamber  tap  was  used  so  the  AOA  information  was  not  believed  to 
be  accurate.  Hence,  no  axial  force  comparisons  with  AOA  are  shown. 
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(C)  Body  with  15“  fins. 


FIGURE  29.  SOME  CONFIGURATIONS  TESTED  AT  M  =  2.01  (REFERENCE  62) 
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HGURE  30A.  STATIC  AERODYNAMICS  FOR  BODY  ALONE  CASE  OF  HGURE  29 
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HGURE  30B.  STATIC  AERODYNAMICS  FOR  10  DEG  FLARE  CASE 
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HGURE  30C.  STATIC  AERODYNAMICS  FOR  15  DEG  FIN  CASE 
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FIGURE  31  A.  WING-BODY-TAIL  CONFIGURATION  CONSIDERED  FOR 
VALIDATION  WITH  AP98  AND  AP95  (REFERENCE  63) 
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HGURE  3  IB.  NORMAL  FORCE  AND  PITCHING  MOMENT  COMPARISONS  OF  THEORY 
AND  EXPERIMENT  FOR  HGURE  31 A  CONHGURATION  (4>  =  0  DEG) 
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The  next  axisymmetric  body  configuration  considered  is  a  22.2  caliber  long  canard-body-tail 
configuration  with  a  100  percent  nose  bluntness  (hemisphere  nose).  The  tail  surfaces  are  fairly  large, 
with  aspect  ratio  0.87,  and  fairly  thick  with  truncated  trailing  edges.  The  canards  are  aspect  ratio  of 
1.73.  The  configuration  is  shown  in  Figure  32A.  The  hangars  which  are  on  the  wind  tunnel  model 
were  not  modeled  by  the  aeroprediction  code.  Tests  were  conducted  for  M„  =  0.2  to  4.63,  AOA  0 
to  20  deg,  control  deflections  0  to  20  deg,  roll  of  0  to  45  deg,  R^/ft  of  2  x  10^  for  a  model  with 
boundary  layer  trips.  Base  pressure  values  as  a  function  of  M„  and  AOA  were  given  in  Reference  64, 
and  these  values  were  added  to  the  axial  force  information  so  total  axial  force  values  could  be 
shown. 


Figure  32B  gives  the  comparison  of  theory  and  experiment  for  ^  =  0  deg  roll  for  both  0  and 
20  deg  control  deflections.  Results  are  shown  in  terms  of  Q,  Cn  and  versus  Mach  number  for 
a  =  0  deg  and  20  deg.  Both  Q  and  are  zero  for  a  =  0  deg  so  only  the  20  deg  control  deflection 
case  is  shown.  Viewing  Figure  32B,  it  is  seen  that  the  AP95  and  AP98  both  give  good  agreement 
to  data  and  are  both  fairly  comparable  in  terms  of  average  accuracy.  The  worst  errors  are  the 
M  =  0.2  pitching  moments  where  theory  gives  a  too  stable  configuration  by  1.6  calibers  or 
7.2  percent  of  the  body  length.  This  is  outside  the  range  of  the  average  accuracy  criteria  of 
±4  percent  of  body  length  for  center  of  pressure.  However,  if  all  Mach  number  values  for  center  of 
pressure  are  averaged,  the  theory  is  well  within  the  ±10  percent  criteria  for  axial  and  normal  force 
and  ±4  percent  for  center  of  pressure. 

The  roll  45  deg  results  are  shown  in  Figure  32C.  Here,  only  experiment  and  AP98  results 
are  given.  Very  good  comparisons  are  given  for  the  AP98  compared  to  data  at  both  control 
deflections  and  both  AOAs  as  a  function  of  Mach  number. 

Two  six-fin  projectile  configurations  are  next  considered  to  show  how  the  AP98  and  AP95 
can  be  applied  to  cases  like  this,  even  though  the  code  only  has  two  and  four  fin  options  at  present. 
The  first  of  these  configurations  is  shown  in  Figure  33A  and  is  taken  from  Reference  65.  It  consists 
of  a  cone-cylinder  body  with  six  tail  fins  for  stability.  The  fins  have  100  percent  blunt  leading  and 
trailing  edges.  Ballistic  range  data  was  available,  and  it  was  assumed  the  model  was  smooth.  Hence, 
the  “wind  tunnel  model  with  no  boundary  layer  trip”  option  was  chosen  in  both  the  AP95  and  AP98 
for  viscous  calculations.  Reynolds  number  was  computed  based  on  sea  level  conditions  at  the  given 
Mach  number  for  a  model  35.2  mm  in  diameter.  Reference  65  gave  both  ballistic  range  and  CFD 
results. 


Figure  33B  shows  the  comparison  of  the  AP95  and  AP98  to  the  ballistic  range  and  CFD 
results  of  Reference  65.  Excellent  agreement  of  the  AP98  and  data  is  shown.  However,  the  AP95 
deviates  from  experiment  significantly  below  M  =  3.25.  The  reason  for  this  large  deviation  is  an 
error  in  the  AP*95  code  for  cases  which  have  large  leading  edge  bluntness.  This  error  was  corrected 
and  is  the  reason  the  AP98  gives  excellent  agreement  with  data  over  the  Mach  number  range  that 
data  was  available.  To  obtain  the  values  from  the  AP95  and  AP98,  the  four  fin  case  was  executed 
with  the  APC.  Then  the  values  in  the  tables  of  data  for  axial  force  of  four  fins  was  multiplied  by  1.5. 
These  values  were  then  added  to  the  body  alone  listed  in  the  tables  of  the  APC.  Hence,  even  though 
the  APC  is  limited  to  four  fins,  six-fin  computations  can  be  obtained  by  using  the  APC  and 
performing  some  hand  calculations  using  the  output  tables  from  the  APC. 
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HGURE  32A.  CANARD-BODY-TAIL  CONHCURATION  WITH  HEMISPHERICAL  NOSE*' 
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TOURE  33A.  FINNED  PROJECTILE  MODEL  (ALL  DIMENSIONS  ARE  IN  CALIBER! 

1  CALIBER  =  35.2  MM) 


HGURE  33B.  AXIAL  FORCE  COMPARISONS  OF  THEORY  TO  BALUSTTC  RANGE 
DATA  FOR  CONHGURATION  OF  FIGURE  33A 
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The  second  of  the  six-jRn  configurations,  and  the  last  axisymmetric  body  case  to  be 
considered,  is  shown  in  Figure  34A  and  the  details  of  the  test  and  CFD  analysis  are  given  in 
Reference  66.  This  configuration  is  quite  similar  to  that  of  Figure  33A  in  that  the  nose  shape  is  an 
8  deg  half  angle  cone  and  the  fin  lea^ng  and  trailing  edges  are  blunt.  The  fins  are  slightly  thicker 
than  those  of  Figure  33A,  however,  and  the  body  is  23.14  versus  13.94  calibers  long  as  well.  The 
diameter  of  the  model  was  27.05  mm  and  the  same  boundary  option  was  chosen  for  the  AP98 
computations  as  for  the  previous  case  of  Reference  65.  Here,  only  AP98  results  are  shown  in 
comparison  to  ballistic  range  and  CFD  results  of  Reference  66. 

Figure  34B  compares  the  axial  force  coefficient,  normal  force  coefficient  derivative  and 
center  of  pressure  predictions  to  ballistic  range  data  and  CFD  results.  Ballistic  range  data  were 
available  for  M„  =  3.5  to  5.3  and  CFD  computations  were  available  at  M  =  5.88, 5.0  and  4.41.  AP98 
results  are  shown  for  M.  =  2  to  6.  Several  points  are  worthy  of  note  here.  First  of  all,  there  is  a 
fairly  high  amount  of  scatter  in  the  ballistic  range  data  for  (Cj^  and  (X(,p)^^Q .  It  is  not  known 
why  this  occurs  but  it  could  be  a  result  of  AOA  motion  that  is  not  accounted  for  in  the  ARC  or  CFD 
computations.  The  AP98  and  CFD  axial  force  results  compare  reasonably  well  with  the  zero  yaw 
drag  results  of  the  ballistic  range  tests  except  at  lower  supersonic  Mach  numbers.  At  M  =  3.5,  the 
ballistic  range  data  shows  a  steeper  increase  in  drag  coefficient  than  does  the  AP98.  It  is  not  clear 
whether  this  is  due  to  an  underprediction  of  the  AP98  compared  to  data  or  yaw  drag  not  accounted 
for  by  the  AP98.  The  previous  six  fin  case  (Figure  33  and  Reference  65)  had  ballistic  range  data  as 
low  as  M  =  3.0.  Excellent  agreement  between  the  AP98  and  data  was  obtained  for  this 
configuration.  Since  the  configuration  is  so  similar,  the  only  plausible  explanation  appears  to  be  yaw 
drag  in  the  data  at  the  lower  Mach  numbers.  C^^  and  X^p  values  of  the  AP98  and  CFD  are  in  very 
close  agreement  and  fall  in  the  middle  of  the  holistic  range  data  given  in  Figure  34B. 


6.2  NONAXISYMMETRIC  BODY  CASES 

The  new  nonaxisymmetric  body  technology  developed  and  described  in  Reference  5  has  been 
evaluated  by  applying  it  to  an  extensive  array  of  aerodynamic  configurations  over  a  broad  range  of 
flight  conditions.  This  included  body  alone  configurations  consisting  of  elliptical,  square  and 
triangular  cross  sections.  The  elliptical  configurations  considered  had  a/b  ratios  varying  from  0.5 
to  2.0.  Square,  diamond  and  triangular  cross  sections  had  comer  roundness  that  varied  between  0 
and  0.33.  Not  all  cases  were  available  on  all  configurations,  however.  Freestream  conditions  varied 
from  Mach  numbers  as  low  as  0.3  to  as  high  as  14.  AOA  as  high  as  60  deg  were  considered. 

Several  wing-body  cases  were  also  considered,^  but  not  nearly  as  much  data  were  available 
as  for  the  body  alone.  Configuration  geometry  and  flight  conditions  were  more  limited  as  well. 
Only  one  wing-body-tail  case^  was  considered  as  this  was  the  only  case  where  data  was  found  in  the 
literature. 
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All  Dimensions  in  Calibers  (One  Caliber  =  27.05  mm) 


HGURE  34A.  SCHEMATIC  OF  M829  PROJECTILE  COISTFIGURATION 
(FROM  REFERENCE  66) 
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HGURE  34B.  COMPARISON  OF  STATIC  AERODYNAMICS  BETWEEN  THEORY 
AND  EXPERIMENT  FOR  HGURE  34A  CONHGURATION 
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Only  three  of  the  Reference  5  results  will  be  shown  here.  All  cases  computed  used  an 
optimum  value  of  the  critical  crossflow  Mach  number  where  transition  from  subcritical  to 
supercritical  conditions  occur.  This  optimum  value  is  not  that  critical  for  higher  Mach  number 
computations  but  is  very  important  to  accurately  predict  subsonic  normal  force.  The  first  case 
considered  is  the  configuration  of  Figure  35.  This  configuration  was  tested  to  24  deg  AOA  at 
=  1.98  and  3.88  in  Reference  50.  All  bodies  in  Figure  35  have  the  same  cross  sectional  area  as 
the  circle.  The  comer  radii  of  the  squares  and  triangles  were  very  small,  so  a  value  of  k  =  0  was 
assumed  in  the  computations.  The  elliptical  shape  10-caliber  body  of  Figure  35  was  tested  iatef*® 
at  Mach  numbers  0.6  to  2.0  and  to  AOA  56  deg.  The  case  shown  here  will  thus  be  the  elliptical  10- 
caliber  long  body  tests  of  Reference  49,  which  go  to  56  deg  AOA,  and  the  square  and  triangular  tests 
of  Reference  50,  which  go  only  to  24  deg  AOA.  Not  all  results  will  be  shown,  as  the  References  50 
and  49  data  bases  were  fairly  extensive.  Most  of  these  results  are  shown  in  Reference  5,  however. 

Figures  36  through  38  give  the  elliptical  body  results  for  Mach  numbers  of  0.6, 1.2  and  2.0. 
Results  shown  are  for  ellipticity  values  of  0.5, 1 .0  and  2.0  and  are  given  in  terms  of  normal  and  axial 
force  coefficients  and  center  of  pressure.  Also,  the  axial  force  coefficient  does  not  include  a  base 
drag  component.  In  examining  the  Figures  36  through  38  comparisons  of  theory  and  experiment, 
it  is  seen  the  theory  does  a  fairly  good  job  of  predicting  the  aerodynamics  on  the  a/b  =  2.0  elliptical 
case  at  all  three  Mach  numbers.  The  normal  force  and  center  of  pressure  predictions  are  quite 
encouraging  as  they  are  well  within  the  average  accuracy  levels  of  ±10  percent  and  ±4  percent  of 
body  length  respectively.  The  axial  force  prediction  comparisons  are  not  as  good  as  desired. 
However,  this  could  be  due  to  measurement  accuracy  where  it  is  difficult  to  get  accurate  axial  force 
measurements  with  a  sting  designed  for  measuring  normal  force  at  high  AOA.  The  a/b  =  0.5  results 
are  not  quite  as  good  as  the  a/b  =  2.0  results,  particularly  at  high  AOA.  This  could  be  due  to  the 
critical  value  of  crossflow  Mach  number  prediction.  The  data  appears  to  support  supercritical  flow 
up  to  AOA  of  56  deg  whereas  the  theory  indicates  the  flow  transitions  to  subcritical  conditions 
around  AOA  36  to  40  deg.  Also  the  center  of  pressure  prediction  for  the  a/b  =  0.5  case  at  low  AOA 
appears  to  indicate  a  center  of  pressure  much  further  toward  the  nose  tip  at  M  =  0.6  than  the  theory 
predicts.  It  is  not  clear  what  mechanism  causes  this.  Since  the  nose  length  is  3  calibers  and  there 
is  no  boattail  present,  intuition  would  lead  one  to  expect  the  center  of  pressure  to  lie  somewhere 
between  1.5  and  3.0  calibers  from  the  nose  tip  near  a  =  0  deg. 

Theoretical  and  experimental  results  for  the  squares  and  triangles  of  Figure  35  are  given  in 
Figures  39  and  40  respectively.  Only  the  10-caliber  long  configuration  results  at  M  =  1.98  are 
shown.  Here  the  results  are  given  in  terms  of  lift  coefficient,  lift  to  drag  ratio  and  center  of  pressure. 
In  general,  comparison  of  theory  and  experiment  for  the  squares  and  diamonds  is  quite  encouraging, 
although  not  as  good  as  the  circular  cross  section  shapes.  The  triangular  shaped  body  predictions 
for  lift  coefficient  tend  to  be  somewhat  low  as  AOA  increases.  However,  this  is  to  be  expected  since 
the  values  for  the  circular  cylinder  results  are  also  low.  Lift  to  drag  ratio  predictions  are  quite  good, 
with  the  peak  values  being  reasonably  well  predicted.  Center  of  pressure  prediction  for  the  triangular 
shape  is  well  within  the  ±4  percent  of  body  length  used  as  a  criteria  for  axisymmetric  bodies. 
However,  the  inverted  triangle  center  of  pressure  predictions  slightly  exceed  this  value. 
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HGURE  35.  BODY  ALONE  CONHGURTIONS^  WITH  ELLIPTICAL,  SQUARE, 
DIAMOND,  TRIANGULAR  AND  INVERTED  TRIANGULAR  SHAPES 
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HGURE  36.  AERODYNAMIC  DATA  FOR  2:1  AND  0.5:1  ELLIPSES  OF  HGURE  35  COMPARED 
TO  CIRCULAR  BODY  AT  M  =  0.6:  (A)  NORMAL  FORCE  COEFFICIENT, 

(B)  AXIAL  FORCE  COEFHCIENT,  (C)  CENTER  OF  PRESSURE 
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HGURE  37.  AERODYNAMIC  DATA  FOR  2:1  AND  0.5:1  ELLIPSES  OF  HGURE  35  COMPARED 
TO  CIRCULAR  BODY  AT  M  =  1.2:  (A)  NORMAL  FORCE  COERFICIENT, 

(B)  AXIAL  FORCE  COEITICIENT,  (C)  CENTER  OF  PRESSURE 
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HGURE  38.  AERODYNAMIC  DATA  FOR  2:1  AND  0.5:1  ELLIPSES  OF  FIGURE  35  COMPARED 
TO  CIRCULAR  BODY  AT  M  =  2.0:  (A)  NORMAL  FORCE  COEFFICIENT, 

(B)  AXIAL  FORCE  COEFFICIENT,  (C)  CENTER  OF  PRESSURE 
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FIGURE  39.  AERODYNAMIC  DATA  OF  SQUARES  (k  =  0.0)  AND  DIAMONDS  (k  =  0.0) 
OF  HGURE  35  COMPARED  TO  CIRCULAR  BODY  AT  M  =  1 .98  («/d  =  10); 

(A)  LIFT  COEFHCIENT,  (B)  LEFT  TO  DRAG  RATIO,  (C)  CENTER  OF  PRESSURE 
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HGURE  40.  AERODYNAMIC  DATA  OF  TRIANGLES  (k  =  0.0)  AND  INVERTED  TRIANGLES 
(k  =  0.0)  OF  HGURE  35  COMPARED  TO  CIRCULAR  BODY  AT  M  =  1 .98  (C/d  =  10):  (A)  LIFT 
COEFHCIENT,  (B)  LEFT  TO  DRAG  RATIO,  (C)  CENTER  OF  PRESSURE 
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The  next  case  considered  is  a  10-caliber  wing-body  and  wing-body-tail  case  shown  in 
Figure  41  (see  Reference  49).  The  body  cross  section  is  an  ellipse  with  an  a/b  =  2.0,  and  the  nose 
length  is  3.0  calibers.  Figures  42  and  43  present  the  normal  force  and  center  of  pressure  comparisons 
of  experiment  and  theory  at  M  =  0.6  and  2.0  for  the  wing-body  and  wing-body-tail  cases 
respectively.  Data  for  the  wing-body-tail  case  at  M  =  2.0  was  only  available  to  a  =  34  deg,  whereas 
all  other  cases  have  data  to  60  deg  AO  A.  As  seen  in  Figures  42  and  43,  theoretical  predictions  are 
quite  good  for  center  of  pressure  and  fair  to  good  for  normal  force  coefficient. 

The  third  configuration  shown  (see  Figure  44)  is  a  waverider  configuration  taken  from 
Reference  67.  Lift,  drag,  and  pitching  moment  comparisons  of  theory  and  experiment  at  M  =  14.0 
to  a  =  25  deg  are  shown  in  Figure  45.  Results  are  quite  encouraging,  even  though  this  configuration 
does  not  quite  fit  the  triangular  shape,  which  has  60  deg  angles  in  all  comers.  Note  that  these  results 
are  based  on  a  375  in.^  planform  area. 


7.0  SUMMARY  AND  FUTURE  OPPORTUNITIES 


In  summary,  new  technology  has  been  developed  to  allow  extension  of  the  aeroprediction 
code  to  the  roll  position  of  45  deg  (fins  in  “x”  or  cross  roll  orientation)  in  addition  to  the  roll  position 
of  0  deg  (fins  in  “+”  or  plus  roll  orientation).  It  has  also  been  extended  to  compute  aerodynamics 
of  nonaxisymmetric  bodies  based  on  an  equivalent  axisymmetric  body.  In  addition,  the  nonlinear 
aerodynamic  loads  have  been  distributed  over  the  body  and  lifting  surfaces  to  provide  a  more  useful 
tool  for  preliminary  structural  analysis.  Finally,  new  technology  was  developed  to  improve  the 
prediction  of  axial  force  at  AOA. 

To  make  the  AP98  more  user  friendly,  an  improved  pre-  and  post-processing,  personal- 
computer  interface  is  being  developed.  These  new  technologies  have  been  integrated  into  the  AP95 
and  will  be  transitioned  to  legitimate  users  as  the  AP98.  Comparisons  of  the  new  theory  have  been 
made  to  both  experimental  data  and  the  AP95.  Comparisons  of  theory  and  experiment  show  the 
AP98  to  be  at  least  as  good  as  the  AP95  and,  in  general,  maybe  slightly  better.  In  general,  average 
accuracy  levels  of  aerodynamics  are  ±10  percent  on  axial  and  normal  force  and  ±4  percent  of  body 
length  on  center  of  pressure.  For  nonaxisymmetric  body  cases,  accuracy  can  be  slightly  higher  than 
these  numbers,  although  not  enough  cases  have  been  considered  to  make  a  definitive  assessment. 
While  these  accuracy  levels  are  encouraging  for  a  semiempirical  code,  they  could  be  improved  upon 
by  using  computational  fluid  dynamics  codes  or  additional  experimental  data  or  both  to  minimize 
assumptions. 

Future  opportunities  are  placed  in  two  categories.  These  categories  are  defined  as  major 
technology  needs  and  value-added  technology  needs.  The  difference  between  the  definition  of 
“major”  and  “value-added”  has  to  do  with  the  cost  to  develop  the  technology  and  integrate  it  into  the 
APC.  Generally,  before  any  new  technology  is  developed,  there  needs  to  be  a  request  from  users  and 
sponsor  support  to  fund  the  technology  development. 
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WING-BODY  WING-BODY-TAIL 

HGURE  41 .  GEOMETRY  OF  THE  WING-BODY  AND  WING-BODY-TAIL  CONHOURATIONS 
WITH  2:1  ELLIPTICAL  BODIES  (FROM  REFERENCE  49) 

60 


FIGURE  42.  AERODYNAMIC  DATA  FOR  THE  WING-BODY  CONFIGURATION  OF  FIGURE  41 
WITH  A  2:1  ELLIPTICAL  CROSS  SECTION  BODY:  (A)  NORMAL  FORCE  COEFFICIENT, 

(B)  CENTER  OF  PRESSURE 
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HGURE  43.  AERODYNAMIC  DATA  FOR  THE  WING-BODY-TAIL  CONFIGURATION  OF 
HGURE  41  WITH  A  2: 1  ELLIPTICAL  CROSS  SECTION  BODY:  (A)  NORMAL  FORCE 
COEFHCIENT,  (B)  CENTER  OF  PRESSURE 


HGURE  44.  WIRE-FRAME  GEOMETRY  OF  THE  WAVERIDER  (FROM  REFERENCE  67) 
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The  NSWCDD  Aeroprediction  Code  has  progressed  to  the  point  that  it  can  compute  planar 
aerodynamics  with  acceptable  accuracy  over  the  configuration  and  flight  envelope  of  interest  to  most 
tactical  weapons.  The  two  remaining  major  technology  needs  are  to  account  for  aerodynamic  effects 
of  side  jets  and  an  accurate  semiempirical  code  for  out  of  plane  aerodynamics.  The  first  problem 
requires  a  good  generic  wind-tunnel  database  and  this  requirement  is  being  marketed  to  BMDO.  The 
latter  need  may  be  beyond  the  scope  of  a  semiempirical  code  due  to  the  complex  out-of-pitch  plane 
aerodynamic  nonlinearities  and  the  fact  that  many  of  these  nonlinearities  are  within  the  accuracy  of 
the  data.  This  is  particularly  tme  for  induced  roll.  As  a  result,  at  NSWCDD,  either  wind  tunnel  data 
or  more  accurate  numerical  codes  have  been  used  to  predict  out  of  plane  aerodynamics.  A  less 
accurate  (±25  percent)  out  of  plane  semiempirical  code  may  be  feasible,  but  it  is  not  clear  if  this 
requirement  is  of  interest  and  general  utility. 

There  are  several  value-added  technology  needs  which  would  enhance  the  capability  of  the 
APC  significantly,  while  being  much  less  costly  than  the  two  major  technology  needs  described 
above.  Some  of  these  include:  a)  using  CFD  to  improve  the  predictions  outside  the  range  of  the  data 
bases;  b)  using  a  new  NASA  data  base  for  fine  tuning  r/s  effects;  c)  integrating  the  APC  with  other 
software  models  to  form  system  engineering  models;  d)  modifying  dynamic  derivative  calculations 
so  flared  vehicles  can  be  handled;  e)  including  a  six-fin  option  into  the  APC;  f)  incorporating  base 
bleed  for  projectile  applications;  g)  including  nonlinear  aerodynamics  into  the  dynamic  derivatives; 
and  h)  incorporating  methodology  for  internal  shock  interactions. 

Over  the  history  of  the  APC  series  (1972  to  the  present),  significant  advances  have  been 
made  in  semiempirical  aeroprediction  methods.  As  seen  in  the  opportunities  list,  much  work  still 
remains  to  be  done. 
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AOA 

APC 

AP95,  AP98 

AR 

DOF 

LT 

MNT 

NASA/LRC 
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SB,  SBT 

Agq 

Ap 

■^REF 

Aw 

a,b 

a' 

b 

Ca 


9.0  SYMBOLS  AND  DEFINITIONS 

Angle  of  Attack 
Aeroprediction  Code 

1995  and  1998  versions  of  the  APC  respectively 
Aspect  Ratio  =  bVAw 
Degrees  of  freedom 
Linear  Theory 

Modified  Newtonian  Theory 

National  Aeronautics  and  Space  Administration/Langley  Research  Center 
Naval  Surface  Warfare  Center,  Dahlgren  Division 
Slender  Body,  Slender-body  Theory 

Cross-sectional  area  of  circular  cylinder  equal  to  that  of  body  with 
noncircurlar  cross  section 

Planform  area  of  the  body  in  the  crossflow  plane  (fF) 

Reference  area  (maximum  cross-sectional  area  of  body,  if  a  body  is 
present,  or  planform  area  of  wing,  if  wing  along)(fP) 

Planform  area  of  wing  in  crossflow  plane  (ft^) 

Semimajor  and  semiminor  axis,  respectively,  of  ellipse 

Body  shape  parameter  (see  Reference  34) 

Wing  span  (not  including  body)(ft) 

Axial  force  coefficient 
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c  c  c 

Ag 


SB 

N 


Base,  skin-friction,  and  wave  components,  respectively,  of  axial  force 
coefficient 

Axial  force  coefficient  at  0  deg  AOA 

Axial  force  coefficient  term  due  to  AOA 

Axial  force  coefficient  term  due  to  control  deflection 

Crossflow  drag  coefficient 

Drag  coefficient 

Circumference  of  body  (ft) 

Lift  coefficient 

Roll  moment  coefficient 

Pitching  moment  coefficient  (based  on  reference  area  and  body  diameter, 
if  body  present,  or  mean  aerodynamic  chord,  if  wing  alone) 

Linear  component  of  pitching  moment  coefficient 

Nonlinear  component  of  pitching  moment  coefficient 

Ratio  of  the  local  normal  force  coefficient  of  a  body  with  a  noncircular 
cross  section  to  that  with  a  circular  cross  section  calculated  by  slender 
body  and  Newtonian  theory  respectively 

Local  normal  force  coefficient 

Side  moment  coefficient 

Normal  force  coefficient 

Normal  force  coefficient  of  body  alone 

Negative  afterbody  normal-force  coefficient  due  to  canard  or  wing-shed 
vortices 

Normal-force  coefficient  on  body  in  presence  of  wing  or  tail 
Linear  component  of  normal-force  coefficient 
Nonlinear  component  of  normal-force  coefficient 
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Cn 

f^T(V) 


Cp 


Cp 

(^PgW.a 

Cy 


Normal  force  coefficient  slope  of  wing  and  tail  respectively 

Negative  normal-force  coefficient  component  on  tail  due  to  wing  or 
canard-shed  vortex 

Normal  force  coefficient  of  wing  alone 
Normal-force  coefficient  of  wing  or  fin  in  presence  of  body 
Normal-force  coefficient  derivative 
Pressure  coefficient  f  — — ] 


l/2p.V: 

\  / 

Base  pressure  coefficient 

Base  pressure  coefficient  of  body  without  fins  at  AOA 
Side  force  coefficient 


Cr 


Rood  chord  (ft) 


Ct 

cal 

dB 

d„f 

^^w(B)  dKg^^ 
da  ’  da 
deg 


Tip  chord  (ft) 

Caliber(s)  (one  body  diameter) 

Body  diameter  (ft)  at  base 
Reference  body  diameter  (ft) 

Rate  at  which  or  ICg^^^  decreases 
Degree(s) 


Pj  C],  C2,  C3 


f,h 


Dimensionless  empirical  factors  used  in  nonlinear  models  of  kw(B)  and 
to  approximate  effects  due  to  high  AOA  or  control  deflection 

Lateral  and  vertical  position  of  wing  vortex 


fw> 


Lateral  location  of  wing  or  tail  vortex  (measured  in  feet  from  body  center 
line) 


i 


Tail  interference  factor 


k,  ki 


Parameters  used  to  define  comer  radius  for  squares  and  triangles 
(k  =  r„AVM;k,=r„AV) 
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Kb(W)>  Kb(T) 

^bCT) 

[^B(W)]mIN 
Kw(B)>  ^T(B) 

kw(B)’  ^^T(B) 

AK 

[AKw(B)]o=o 

and 

a=0 

^ref 

m 


Mn 

M 


Nr 


M. 


M. 

NF 

NFo 

NF, 

P 

Pw 

r 


Ratio  of  additional  body  normal-force  coefficient  in  presence  of  wing,  or 
tail  to  wing,  or  tail  alone  normal-force  coefficient  at  6  =  0  deg 

Ratio  of  additional  body  normal-force  coefficient  due  to  presence  of  wing 
or  tail  at  a  control  deflection  to  that  of  wing  or  tail  alone  at  a  =  0  deg 

Minimum  value  of  Kb(W)  as  percent  of  slender-body  theory  value 

Ratio  of  normal-force  coefficient  of  wing  or  tail  in  presence  of  body  to  that 
of  wing  or  tail  alone  at  6  =  0  deg 

Ratio  of  wing  or  tail  normal-force  coefficient  in  presence  of  body  due  to  a 
control  deflection  to  that  of  wing  or  tail  alone  at  a  =  0  deg 

Nonlinear  component  of  wing-body  or  body-wing  interference 

Amount  that  the  experimental  values  of  K^cb)  and  Kb(W)  exceed  slender 
body  theory  at  a  =  0  deg 


Reference  length  (ft) 

Cotangent  of  leading  edge  sweep  angle 

Mach  number  normal  to  body  =  M„  sin  a 

Normal  Mach  number  where  flow  transitions  from  subcritical  to 
supercritical  conditions 

Value  of  normal  Mach  number  to  body  used  in  Newtonian  correction 
factor  calculation 

Freestream  Mach  number 

Newtonian  correction  factor 

Newtonian  correction  factor  for  crossflow  Mach  number  of  zero  on 
squares  and  triangles 

Newtonian  correction  factor  for  an  ellipse  at  a  ^  20  deg 
Pressure  (Ib/ft^) 

Loading  factors  in  leeward  and  windward  planes  respectively 
Local  body  radius  (ft) 
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Tgqj  d 


eq 


R. 


s 


SBTSF 

V„ 


''N 

w 

w„ 


Xcp 


(^cp)l>  (^cp)i 


'ML 


x,y,z 

a 

OCc 


ttp 


“m 


Radius  and  diameter,  respectively,  of  a  circular  cross-section  body  which 
has  same  cross-sectional  area  as  that  of  noncircular  cross-section  body 

Comer  radius  of  a  rounded  comer  on  square  or  triangle 

Re5molds  number  where  flow  transitions  from  subcritical  to  supercritical 
conditions 

Re5molds  number  based  on  body  diameter 

An  effective  Reynolds  number  above  which  the  flow  transitions  from 
subcritical  to  supercritical  conditions 

Radius  of  body  at  wing  or  tail  locations 

Wing  or  tail  semispan  plus  the  body  radius  in  wing-body  lift  methodology 

Slender  body  theory  scaling  factor 

Freestream  velocity 

Velocity  normal  to  body 

Length  of  one  side  of  a  triangle  or  square 

Maximum  diameter  of  a  triangle  or  square  as  measured  normal  to  the 
velocity  vector 

Center  of  pressure  (in  feet  or  calibers  from  some  reference  point  that  can 
be  specified)  in  x  direction 

Center  of  pressure  of  linear  and  nonlinear  terms  of  normal  force 
Axis  system  fixed  with  x  along  centerline  of  body 
Angle  of  attack  (deg) 

Angle  of  attack  where  wing-body  interference  factor  starts  decreasing 
(deg) 

Angle  of  attack  where  the  wing-body  interference  factor  reaches  a 
minimum  (deg) 

Angle  of  attack  where  reaches  a  constant  value 
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Local  angle  of  attack  of  wing  or  tail  (a^  +  6  or  %  +  respectively,  in 
degrees) 

Angles  of  attack  used  in  nonlinear  model  for  Kgjw) 

-  1 

Control  deflection  (deg),  positive  leading  edge  up 

Angle  between  velocity  vector  and  tangent  to  body  surface 

Deflection  of  wing  or  tail  surfaces  (deg),  positive  leading  edge  up 

Parameter  used  in  viscous  crossflow  theory  for  nonlinear  body  normal 
force  (in  this  context,  it  is  the  normal  force  of  a  circular  cylinder  of  given 
length-to-diameter  ratio  to  that  of  a  cylinder  of  infinite  length) 

Decay  constant  used  in  second-order-shock-expansion  theory 

Vortex  circulation,  positive  counterclockwise  facing  upstream  (ftVsec) 

Roll  position  of  missile  fins  (0  =  0  deg  corresponds  to  fins  in  the  plus  (+) 
orientation).  O  =  45  deg  corresponds  to  fins  rolled  to  the  cross  (x) 
orientation 

Taper  ratio  of  a  lifting  surface  =  c/c, 

Roll  position  of  point  on  body  with  <J)  =  0  deg  being  the  leeward  plane 
Local  surface  slope  of  body  with  respect  to  body  axis 


Freestream  conditions 
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